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Chapter  1 
Introduction 


There  is,  as  shown  in  Fig.  1.1,  a  circular  waveguide  which  is  closed  at  one  end. 
Two  symmetrically  placed  apertures  in  the  lateral  wall  of  this  waveguide  are 
backed  by  rectangular  waveguides  of  identical  dimensions.  The  dimensions  of 
the  waveguides  are  such  that  only  the  TEU  and  TA/01  modes  can  propagate 
in  the  circular  waveguide  and  that  only  the  TE\o  mode  can  propagate  in 
the  rectangular  waveguides.  The  problem  is,  as  stated  in  [1,  Chapter  1],  to 
find  out  how  much  of  the  power  of  an  incident  T M0 j  wave  in  the  circular 
waveguide  is  reflected  in  the  circular  waveguide  and  how  much  of  this  power 
is  transmitted  into  the  rectangular  waveguides. 

In  this  report,  the  analytical  results  of  [1]  for  this  problem  are  manipulated 
into  expressions  suitable  for  evaluation  by  means  of  a  digital  computer.  These 
analytical  results  are  not  derived  here;  they  are  merely  referred  to.  For  this 
reason,  the  reader  of  this  report  should  obtain  a  copy  of  reference  [I]. 

A  computer  program  was  written  in  FORTRAN.  Some  numerical  re¬ 
sults  obtained  by  using  this  computer  program  are  presented.  The  computer 
program  will  be  described  and  listed  in  a  forthcoming  report. 

As  is  shown  in  Fig.  1.2,  the  interiors  of  the  left-hand  rectangular  wave¬ 
guide,  the  right-hand  rectangular  waveguide,  and  the  circular  waveguide  are 
called  regions  1,  2,  and  3,  respectively.  The  electromagnetic  field  in  region  1  is 
radiated  by  The  field  in  region  2  is  radiated  by  M}2'1-  The  field  in  region 
3  is  radiated  by  the  combination  of  ,  —A/*1*.  and  —  A/^2'.  The  magnetic 
currents  in  Fig.  1.2  are  supposed  to  be  loacated  right  on  (infinitesimal  dis¬ 
tances  from  either  side)  of  the  closing  conductors.  The  finite  displacement 
of  these  magnetic  currents  from  the  closing  conductors  in  Fig.  1.2  is  only  for 


1 


REGION  1 


z 


X 


REGION  2 


Fig.  1.1.  Top  and  side  views  of  the  T A/01  to  T El0  mode  converter 


the  purpose  of  illustration.  The  magnetic  currents  A/*1*  and  A/(2*  are  given 
by  [1,  eqs.  (2.11)  and  (2.12)]  in  which  the  l/’s  are  the  elements  of  the  column 
vector  on  the  left-hand  side  of  [1,  eq.  (2.22)].  In  [1,  eq.  (2.22)],  V1 .  V'2,  and 
V'3  are  the  admittance  matrices  for  regions  1,  2,  and  3,  respectively.  The 
column  vector  on  the  right-hand  side  of  [l,  eq.  (2.22)]  is  called  the  excitation 
vector.  The  matrices  Yx  and  V-2  are  treated  in  Chapter  2,  the  matrix  Y3  is 
treated  in  Chapter  3,  and  the  excitation  vector  is  treated  in  Chapter  4. 

After  solving  [1,  eq.  (2.22)]  for  the  E’s  which  determine  A/(1*  and  A/^ 
according  to  [l,  eqs.  (2.11)  and  (2.12)],  we  find  the  electromagnetic  fields  in 
regions  1,  2,  and  3.  Expressions  for  the  fields  in  regions  1  and  2  are  obtained 
in  Chapter  5.  The  field  in  region  3  due  to  the  combination  of  —A l}1'1  and 
—  A^2)  is  considered  in  Chapter  6.  An  expression  is  obtained  for  this  field 
below  the  apertures  where  z  <  —c/2.  Expressions  are  also  obtained  for  this 
field  in  the  apertures.  In  Chapter  7,  numerical  results  are  presented  for  the 
power  transmitted  into  the  rectangular  waveguides,  the  power  reflected  back 
into  the  circular  waveguide,  and  the  magnitudes  of  the  6-  and  2-components 
of  the  electric  field  in  one  of  the  apertures  when  a  T Eq\  wave  is  incident 
in  the  circular  waveguide  and  when  the  loads  Z\  and  Z2  that  terminate  the 
rectangular  waveguides  are  both  matched  loads. 

In  Appendix  A,  the  expansion  functions  { A/ -  M.2Jn  } 
are  ordered  so  that  each  one  of  them  can  be  identified  by  means  of  a  sin¬ 
gle  positive  integer  rather  than  the  attached  combination  of  subscripts  and 
superscripts.  In  Appendix  B,  a  numerical  procedure  for  obtaining  roots  of 
Bessel  functions  and  their  derivatives  is  described.  Heretofore,  the  excita¬ 
tion  has  simply  been  a  2-traveling  TA/qj  wave  in  the  circular  waveguide.  In 
Appendix  C,  the  response  due  to  this  excitation  is  used  to  find  the  response 
due  to  excitation  by  a  transverse  sheet  of  TMqi  electric  current  between  two 
impedance  loads  as  shown  in  Fig.  C.l. 
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Chapter  2 

The  Admittance  Matrices  for 
the  Rectangular  Waveguides 


The  admittance  matrice  for  region  1,  the  left-hand  rectangular  waveguide  in 
Fig.  1.2,  is  V'1  given  by  [1,  eq.  (2.25)]  where  the  Y's  on  the  right-hand  side  of 
[1,  eq.  (2.25)]  are  approximated  by  V”s  given  by  [1,  eqs.  (3.44)-(3.47)].  The 
superscripts  on  the  V’s  are  the  same  as  those  on  the  Y's. 

For  convenience,  [1,  eqs.  (3.45)  and  (3.46)]  are  repeated: 

yl.lTE.lTM  _  q 

Y>''™'lTE  =  Q 

In  [1,  eqs.  (3.44)  and  (3.47)],  we  have  [1,  eq.  i 

{j  flpq  l  kp? 

K, 

where  k  =  in  which  ui  is  the  angular  frequency.  Moreover,  fi  and  t  are. 

respectively,  the  permeability  and  permittivity  of  the  medium  in  regions  1. 
2.  and  3.  In  (2.3), 


;a.i2)] 

<  k 
>  k 


(2.1) 

(2.2) 


(2.3) 
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(2.4) 

(2.5) 


There  is  a  correspondence  between  each  pair  of  integers  (p,  q)  used  in  [1, 
eqs.  (3.44)  and  (3.47)]  and  the  subscript  j  in  [1,  eqs.  (3.44)  and  (3.47)]. 
This  correspondence  is  described  in  Appendix  A.  If  (p,  q)  =  (1,0),  then, 
because  the  T E\o  mode  propagates  in  the  rectangular  waveguide,  k  >  fc10 
and  substitution  of  (2.3)  into  [1,  eq.  (3.47)]  and  subsequent  multiplication  by 
—jrj  where  p  =  gives 


-  *>1 ATEATE 

-JvYij 


■  Ao(coq  floSi  +  jZxYfif  sm  fax  j) 
k(j  sin  flioxi  +  Z\YiqE  cos  /310^i)  U' 

(P,9)  =  (l,0)  (2.6) 


where 


x  j  —  L\  x  q 


a  sin  (f)Q 


<f>0 


(2.7) 

(2.8) 

(2.9) 


Yx T0E  is  the  characteristic  admittance  of  the  T E\q 
waveguide  [l,  eq.  (A. 25)] 

VTE  _  A o 

Y'°  “  hr,' 

and  8{j  is  the  Kronecker  delta  function  given  by 


mode  in  the  rectangular 


(2.10) 


8  ij 


1 ,  i  =  j 
i^j 


(2.11) 


The  identities  [2,  formulas  65^  6  and  654.7]  were  used  to  obtain  (2.6).  In 
(2.6),  the  subscript  j  is  not  to  be  confused  with  the  other  j's.  Each  of 
these  other  j's  is  \J—\.  Because  the  T Epq  mode  does  not  propagate  in  he 
rectangular  waveguide  when  (p,  q)  /  (1,0),  ~fpq  is  real  when  (p,  q)  ^  (l,o) 
and,  from  [l,  eq.  (3.47)], 

<P.»)  +  (1-0).  (2.12) 

The  factor  —  j  on  the  left-hand  side  of  (2.12)  has  rendered  the  right-hand 
side  of  (2.12)  real.  The  factor  p  on  the  left-hand  sides  of  (2.6)  and  (2.12)  has 
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rendered  the  right-hand  sides  of  (2.6)  and  (2.12)  independent  of  rj.  Multi¬ 
plication  of  [1,  eq.  (3.44)]  by  —jr)  gives 


_  J^c 


(2.13) 


Because  all  the  TM  modes  in  the  ctangular  waveguide  are  evanescent,  7P, 
is  real  in  (2.13)  so  that  the  right-hand  side  of  (2.13)  is  real. 

The  admittance  matrix  for  region  2,  the  right-hand  rectangular  waveguide 
in  Fig.  1.2,  is  Y 2  given  by  [1,  eq.  (2.27)]  where  the  Y's  on  the  right-hand  side 
are  approximated  by  T’s  given  by  [1,  eqs.  (3.49)— (3 .52)] .  Similar  to  (2.1), 
(2.2),  (2.6),  (2.7),  (2.12),  and  (2.13),  we  have 


y2,2TE,2TM 

y2,2TM,2TE 

x'i 

.jr)Y^TE,2TE 


=  0 
=  0 

_  ,/9io(cos  0lox2  +  jZ2Y?0E  sm0lox2)  ^ 

J  k(j  sin  (3l0x2  +  Z2Y&e  cos  (310x2 )  '3' 

(p-i  q)  =  (i,o) 


x2  =  l2  -  x0 


■  Ss2,2TM,2TM 

-JViYij 


{p,q)  ±  (1,0) 


(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 

(2.19) 
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Chapter  3 


The  Admittance  Matrix  for 
the  Circular  Waveguide 


The  admittance  matrix  for  region  3,  the  circular  waveguide  in  Fig.  1.2,  is  Y3 
given  by  [1,  eq.  (2.29)]  where  the  Y's  on  the  right-hand  side  of  [1,  eq.  (2.29)] 
are  given  by  [1,  eqs.  (4.93),  (4.111),  (4.112),  and  (4.113)]  in  which  T,  Su  S2, 
53,  5j,  and  55  are  given  by  [1,  eqs.  (4.94)— (4.99)].  The  previously  mentioned 
equations  for  the  T’s  are  recast  as 

oo  oo 

- JnY, I*™"™  =  £  £  T{W85i  +  W9S3  -  Wl0S<  -  WnSs}  (3.1) 

r=0  3=1 
00  00 

- M H*™*™  =  ££ f{vv9s, - - WuS,  +  vv.oS,}  (3.2) 

r= 0  3=1 

oo  oo 

-jriY^a™'yTE  =  EEH  &  +  WnS3  +  WsS<  +  \ F953}  (3.3) 

r= 0  3=1 

oo  oo 

-jvY,3J'aTE'yTE  =  'EHfWiSx  -  VF1053  +  W9s4  -  1F855}  (3.4) 

r=0  3=1 

where 


Wio 

Wll 

.S\ 


5,  = 


54 

55 

^72 

^73 

^7i 


Z\  = 

r2  = 

-3  = 

Z4  = 

‘5  — 

;(i)  _ 

,(2) 

i(3) 

-('!) 

4(5) 


np 
mpc 

~r 

(Zl  -  Z2 )<^2 

z4^3 

Zs 4>^ 

(_  l)*^1 

^72 

(_1)T-^3 
(-l)a+7<^a74 

MW^0V;U™«)x(1)  /er' 


XraJr+l(Xrs) 

r2Jl{klEa)zM 
(X'r]  -  r2)J2(4j 
{kJsEa)2rJ?(kJEa)zW  /sin  0O 


(Xr«  -  r2)Jr2«s)  V  <j>0 

( kJsEa)2rJr(kJ,Ea)z{4)  ( sin  0O 
(x;2-r2)J2(*;j 

(kJsEa)4Jr(kJ3Ea)z{5)  ( sin  4>, 

(Xr«  ~  r2)Jr{X'rs)  \ 

4z<!) 


c2-<;JsXI  a 

4^S*. 


4-<3) 


c2 

4-(5) 

c27rr3E< 


(3-8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 

(3.22) 

(3.23) 

(3.24) 

(3.25) 

(3.26) 

(3.27) 
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In  (3.5)-(3.9),  j  determines  p  and  q  in  the  manner  described  in  Appendix 
A;  i  determines  m  and  n  in  the  same  manner.  In  obtaining  (3.1)-(3.4),  we 
started  the  index  r  of  summation  in  [1,  eqs.  (4.96)  and  (4.97)]  at  0  instead 
of  1.  This  was  possible  because  the  r  =  0  terms  so  introduced  are  zero.  We 
multiplied  [1,  eq.  (4.93)]  by  <Jtmentpeq/ 16,  [1,  eq.  (4.111)]  by  yfe^ej 4,  and 

[1,  eq.  (4.112)]  by  \JemtnlA.  This  was  possible  because  none  of  the  indices 
m,  n,  p,  and  q  which  appear  in  these  multipliers  is  ever  zero. 


3.1  Evaluation  of  the  0’s  in  (3.10)— (3.13) 

The  6's  in  (3. 10)— (3. 13)  are  given  by  (3. 14)— (3 . 17) .  Substitution  of  [1, 


eqs.  (E.31  )-(E.34)]  into  (3. 14)— (3. 1 7)  gives 

(3. 28) 

<^2  =  42Va1^2  +  ^pl)4>alyl  (3.29) 

<^3  =  (3.30) 

^  =  (-1)®+^  ^(4)^2  +  .  (3.31) 

For  7  =  a,  we  have  [1,  eqs.  (E.46)-(E.49)] 

P"'  =  >,  7  =  a  (3.32) 

P™  =  *£>,  7  =  o  (3.33) 

P™  =  7  =  0  (3.34) 

P ™  =  *£’,  7  =  0-  (3.35) 

For  7  ^  a,  we  have  [1,  eqs.  (E.53)-(E.56)] 

_  (_i)r  j<^>sin  —  -  <j>^cos  —  1 ,  7/a  (3.36) 

4>a2^  =  (_l)r  j</£>sin—  ~<t>{*]cos~  1,  7/a  (3.37) 

l  x0  j 

<£“H2  =  (-l)r  {^cos  —  +  ^sin— ]  ,  7^0  (3.38) 

l  x0  x0) 

4>a2~<2  =  (-l)r  j<^)Cos  —  +  ^g)sin  —  1,  7 /a.  (3.39) 

{  x„  x0J 
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(3.40) 


In  (3.2S)— (3.31),  the  <f>p's  are  given  by  [1,  eqs.  (E.23)-(E.26)] 


where 


^2)  = 

= 


b  (  sin  A  sin  A+  ] 

2^1“^  A+~  J 

b  fsin2(A“/2)  sin2(A+/2)l 

2^  {  (A- 12)  +  (A+/2)  J 

b  f  sin2(A~/2f  sin2(A+/2) ) 
toll  (.4-/2)  +  (A+/2)  j 

b  (  sin  A~  sin  A+  ] 

2^1  A-  +  A+  J 


..  rb 

A  =  pit  H - 

Xq 


(3.41) 

(3.42) 

(3.43) 

(3.44) 

(3.45) 


The  0’s  in  (3.36)— (3.39 )  are  given  by  the  right-hand  sides  of  (3.40)~(3.43) 
with  A+  and  A-  replaced  by  the  right-hand  sides  of  (3.44)  and  (3.45)  with 
p  replaced  by  m. 


3.2  Evaluation  of  the  z*  s  in  (3.10)— (3.13) 

The  z’s  in  (3.10)— (3.13)  are  given  by  (3.1S)— (3.22).  In  (3.1S)— (3.22),  er  is 
Neumann’s  number  given  by  [1,  eq.  (B.9)j 


1,  r  =  0 

2  r  =  1  2  •  •  •  ’ 

“7  '  L  7  “7 


(3.46) 


JT  is  the  Bessel  function  of  the  first  kind  of  order  r,  xTa  is  the  root  of  Jr, 
J'r  is  the  derivative  of  JT  with  respect  to  its  argument,  and  x'TS  is  the  root 
of  J'.  The  roots  (xrj,  r  =  s  =  1, 2,  •  •  •}  and  {x'T3,  r  =  1, 2,  ■  •  •}  are  ordered 
such  that 


0  <  Xrl  <  xr2  <  Xr3  •  •  •  (3.47) 

0  <  x'rX  <  xr2  <  x'r3  •  •  •  (3.48) 

II 


Still  in  (3. IS)— (3.22),  we  have  [1,  eqs.  (B.7)  and  (B.41)] 

.  TM  _  £rs 

Krs  ~  „ 

a 

l.TE  _  ^rs_ 

Kr s  ~  „  ■ 

a 

Since  Jr{xrs)  =  0,  we  have  [3,  formula  9.1.27] 

Jr(xTs )  Jr+l{xrs)- 


(3.49) 

(3.50) 


(3.51) 


Substitution  of  (3.49)— (3.51 )  into  (3.1S)-(3.22)  yields 


r2z(2) 


_  rx^i(3)  / sin  ^0\ 
~3  x'rl  -  r2  \  4>0  ) 


(3.52) 

(3.53) 

(3.54) 

(3.55) 

(3.56) 


The  z’s  in  (3.52)-(3.56)  are  given  by  (3.23)-(3.27)  where  [1,  eqs.  (B.24) 
and  (B.53)] 


where 


>  xn  <  ka 

yjx2,  -  ( ka )2  ,  xr,  >  ka 
]3j3Ea  ,  x'rs  <  ka 

y/K3.  ~  CkaY  ,  4,  >  ka 


=  \](ka)2  -  x2TS 
Pj,Ea  =  \J(ka)2  -  x£. 


(3.57) 

(3.58) 
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(3.59) 

(3.60) 


From  [1,  eqs.  (F.60),  (F.61),  (F.76),  (F.81),  (F.118),  and  (F.119)], 


-d)  =  l  (DTMGTM  +  C2FTM\ 

4  v  / 

(3.61) 

-(2)  _  l  (^DTEqTE  +  JfTE'j 

(3.62) 

.(3)  =  1  (D(3)G(3,  +,c2F(3)^ 

(3.63) 

_  _I  ( DteG(4)  +  C2F(4)) 

(3.64) 

TE  TE 

_(5)  _  Irs  Zss  |  •?  / F(3)G(4)  ,  C2F(5)\ 

(kjf)2  4  V  ) 

(3.65) 

where  the  D' s,  the  G”s,  the  F’s,  and  zTE  are  dealt  with  in 

[1,  Appendix  F]. 

In  view  of  [l,  eq.  ( F .83)]  and  (3.50),  substitution  of  (3.61 )— (3 .65)  into  (3.23)- 

(3.27)  gives 

i  (F™  +  6™r>™) 

(3.66) 

:m  =  j  (f  TE  +  dTEb TE)  T«0 

(3.67) 

;(3|  _  f'lJ)  +  gTE Dm 

(3.68) 

i'4>  =  -  (F«  +  C<4’£>™) 

(3.69) 

-(si  4«4£  ,  j  (f(5)  +  gW) 

(3.70) 

where 

Gf  =  — ,  8  =  TM,  TE 
q  c 

(3.71) 

ns 

Dsn  =  — ,  8  =  T M,  TE 

c 

(3.72) 

*. - 0131 
c 

(3.73) 

GG) 

Gf  =  — . 

1  c 

(3.74) 

13 


3.2.1  The  TM  Quantity  for  xrs  <  ka 

With  xrs  <  ka,  substitution  of  (3.57)  into  (3.66)  gives 


M1)  =  F™  +  GV'b™ 

Substituting  [1,  eq.  (F.37)]  into  (3.71),  we  obtain 

sin{q5-c)cos(f3srsL$)  -  2  sin2  sin(^slj ) 


(3.75) 


3  = 


qs~c 


sin(r/5+c)  cos(/3^Z,J )  4-  2  sin2 


qs+cs 


sin  (0r.Lt) 


+• 


qs+c 


(3.76) 


where  8  may  be  either  TM  or  TE  and  where  [1,  eqs.  (F.ll),  (F.23),  and 
(F -24)] 


Lt  =  L,  +  C- 


J+.  _ 


q  C  =  qx  -  /3°ac 

3S 

Jrs' 


q  +c  =  qx  +  f3°,C . 


(3.77) 

(3.78) 

(3.79) 


Using  (3.77)— (3.79)  and  [2,  formulas  403.02,  401.03,  and  401.04],  we  reduce 
(3.76)  to 


sin 


yv 


sin 


'q5~c' 


qs+c 


+  (-l)q 


q6~c 


>  cos 


(A*--? 


(3.80) 


2  2 

Substituting  [1,  eq.  (F.S5)]  into  (3.72),  we  obtain 


D*n  = 


j -js\n{ns~c)  -  2 sin2  (!iy£)  —j  sin (n5+c)  +  2  sin2  (—s) 


ns  c 


-itt.Lt 


ns+c 


(3.81) 
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where  nd  c  and  ns+  c  are  given  by  (3.78)  and  (3.79)  with  q  replaced  by  n.  In 
the  same  manner  as  we  reduced  (3.76)  to  (3. SO),  we  can  reduce  (3.81)  to 


Di 


As  for  the  quantity  FT,U  in  (3.75),  we  have  [1,  eq.  (F.79)] 

F8  =  —f(ns~c,  —qS~c)  +  f(ns+c,  qS~c) 
-  f{ns~c ,  qs+c)  +  f(n5+c,  - qS+c ) 


(3.82) 


(3.83) 


where  8  may  be  either  TM  or  TE  and  [1,  eq.  (F.97)] 

sin  a:  j  x  +  y  /  0 

~  y*  '  \  bl  >  f 

(-l)'sin  y  j  x+y  ^  0 

Isl  <  f 

x  +  y  =  0 

|y|  >  0-i 

y_  _t_  (  x  +  y  =  0 

.  3!  5!  +  7!  ’  \  |y|  <  0.1 

where  I  is  the  integer  that  satisfies 


/(•F2/)  =  < 


yx 

y  -  sin  y 
V 2 


’  l 


(3.84) 


x  +  y  =  In .  (3.85) 

When  xr3  <  ka ,  the  TM  quantity  is  now  given  by  (3.75)  in  which 
fiJXIa,  GjX{,  D™ ,  and  F™  are  given  by  (3.59),  (3.80),  (3.S2)  and  (3.83), 
respectively. 


3.2.2  The  TM  Quantity  for  xrs  >  ka 

When  xr3  >  ka ,  we  proceed  to  evaluate  expression  (3.66)  for  r*1),  which  con¬ 
tains  the  quantities  G^XI .D™ ,  and  FTXI .  These  quantities  with  the  super¬ 
script  TM  replaced  by  8  are,  according  to  (3.71),  (3.72),  and  [1,  eqs.  (F.104), 
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(F.100),  (F.101),  and  ( F. l07)-( F.  109 )] ,  given  by 


where 


Gt  = 


<? 

Dt 


0  j c  t-  3(^3  ~  j ) 


=  ^ jcne 


F=j< 


Cq  {sinh(^X+)  -  (  — l)9  sinh  -  c))} 

—  sinh 
cosh 

f  (  — l)ncnc,sinh(<7c) 

(-l)ncnc,  sinh(<7c)  -  cg  , 

4{sinh(<?c)  -  gc] 


[gcf 


n  even 
n  odd 
q^n 
q  =  n  0 

q  =  n  =  0 


9  =  7r 


9-/ 


cn  — 


(  n  7T  )2  +  (7r5C)2 


(3.86) 

(3.87) 


(3.88) 


(3.89) 

(3.90) 


Furthermore,  c,  is  the  right-hand  side  of  (3.90)  with  n  replaced  by  q.  The 
truncated  series  approximation  inherent  in  [1,  eq.  (F.110)]  is  introduced  later 
in  this  section.  The  case  where  iTS  =  ka  is  not  allowed  because,  if  xT,a  =  ka, 
then,  according  to  (3.57),  7 J*ra  would  be  zero  so  that  division  by  7 ™a  in 
(3.66)  would  be  impossible. 

Substituting  (3.77)  into  (3.86)  and  using  [2,  formulas  651.06  and  651.07], 
we  obtain 

M  _  o  /  sinhfoc/2)  cosh(^L3)  ,  q  even 
q  ~Cq  (  sinh(</L3)  cosh(gc/2)  ,  q  odd 

Substitution  of  (3.77)  into  (3.87)  gives 

—  sinh(<7c/2)  ,  n  even  1  H  991 

cosh(<7c/2)  ,  n  odd  J 

Combining  (3.SS),  (3.91),  and  (3.92)  and  using  [2,  formulas  652.12  and  654.5] 
to  simplify  the  result,  we  obtain 


D‘n  =  2jcne-^>  | 


(3.91) 


j(F‘  +  Gl1bi)  = 


~nq 


4"  Cn 


1 

0 


7  =  n  ^  0  | 

otherwise  j 


(3.93) 
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where 


Cfiq  — 


CnCq~ee  ? 


CnCqZee  , 


CnCqZoe  ? 


CnCq^oe  i 


CnCq~oo 


■{: 


n  even 
q  even 
q^n 
n  even 
q  =  n  ^  0 
q  =  n  =  0 
n  odd 
q  even 

n  even 
q  odd 

n  odd 
q  odd 


(3.94) 


In  (3.94), 


=  2  je  29^3  sink  ~  —  e  ^  j  sink  ~ 


~o  —  -ee  4”  9*~ 
~  ~  P-I9L3 


coe  =  —  e  33  sinh  gc 

z00  —  2  je-291,3  cosh  ~  —  e-^  j  cosh  —■  . 


(3.95) 

(3.96) 

(3.97) 

(3.98) 


The  value  of  sinh  {gc/2)  is  excessively  large  when  gcj 2  is  only  moderately 
large.  Moreover,  z0  approaches  zero  more  rapidly  than  gc  as  gc  approaches 
zero.  To  avoid  computational  difficulties,  we  replace  (3.95)-(3.98)  by 


e~9C  _  e-Z9^3  _  1  +  i  |e-3(2L3~c)  +  e-3(2Z,3+c)|  ,  gC  >  l 


2  {e-29^3  sinh  ^  j  sinh  ^  ,  #c  <  1 

e-3c  _  e-2fl^3  -1  +  4  |e-3(2/.3-c)  +  e-3(2i,3+c)|  +  gc ,  (JC  >  J 

2  je-29^3  sinh  ^  —  e-  V  j  sinh  ^  +  gc  ,  0.01  < 

2^  (2e~^  sinh  ^  +  e-29^3  sinh  sinh  ^ 


(3.99) 


0.01  <  gc  <  1 


(3.100) 


((?)  \(?) .(¥)■)  \ 
■("3r  +  ~  +  “7rj  j 


,  gc  <  0.01 
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±  |e-3(2i-3+c)  _  e-3(2L3-c)|  ,  gc>  1 

— e~2aL3  sinh  gc  ,  gc  <  1 

e"2^3  _  e-gc  -  1  +  I  |e-3(2L3-c)  +  e-g(2L3+c)  >  j 

2  je-2^3  cosh  ^  —  e-^}  cosh  ^  ,  #c  <  1 

We  used  the  approximation  [2,  formula  657.1] 

x3  x°  x‘ 

sinh  x  —  x=—  +  —  +  — 

3!  o!  (! 


(3.101) 

(3.102) 


(3.103) 


to  obtain  zQ  of  (3.100)  for  gc  <  0.01.  From  Fig.  1.2.  L 3  >  c/2  so  that  all 
the  arguments  of  the  exponentials  in  (3.99)— (3.102)  are  less  than  or  equal 
to  zero.  The  exponentials  in  ( 3.99 )— (3. 102 )  will  be  excessively  small  when 
gL%  is  only  moderately  large.  However,  this  will  not  cause  any  difficulty  if 
we  use  a  computing  system  which  treats  an  underflow  by  setting  the  number 
equal  to  zero  and  proceeding  without  an  error  message.  The  right-hand 
side  of  (3.94)  cannot  be  evaluated  when  gc  =  0  and  when  n  =  0  or  q  =  0 
because  cn  =  2 /(gc)  when  n  =  0  and  cq  =  2 /(gc)  when  7  =  0.  However,  the 
right-hand  side  of  (3.94)  remains  finite  as  gc  approaches  zero  when  n  =  0  or 
<7  =  0. 

When  xrs  >  ka,  the  TM  quantity  £,l)  is  now  given  by  (3.66)  in  which 
7 J3fa  and  j(FT3r  +  Grq31  D*sl)  are  given  by  (3.57)  and  (3.93),  respectively. 
In  (3.93),  cn  and  cnq  are  given  by  (3.90)  and  (3.94),  respectively.  In  (3.94), 
zee,  ~0,  zoe ,  and  zoa  are  given  by  (3.99)-(3. 102)  where  g  is  given  by  (3. 89). 


3.2.3  The  TE  Quantities  £(2)_£(;>)  for  x'rs  <  ka 

When  x'T3  <  ka,  7 J3Ea  is  given  by  (3. 58)  so  that  expressions  (3.67)-(3.70)  for 
£(2)_~(o)  become 


,<2>  =  _  (f  TE  +  qTE  jjTE'j  3TEa 

z( 3)  =  F(3)  4-  GTqED{3) 
i(4)  =  -  (F(4)  +  G[A)DTnE) 


:(5)  _ 


F(5)  +  G<4)£)t 3) 

0lEa 


+ 


sa 


cx 


<2 


n  =  <7  7^  0 
otherwise 


(3.104) 

(3.105) 

(3.106) 


(3.107) 


IS 


0 


We  obtained  (3.107)  by  substituting  [1,  eq.  IF. 121)]  for  zJE. 

In  (3.104),  fiJEa,  GEE ,  D„E,  and  FTE  are  given  by  (3. GO),  (3. SO),  (3.82), 
and  (3.83),  respectively.  As  for  G^  in  (3.106)  and  (3.107),  substitution  of 
[1,  eq.  (F.  122)]  into  (3.74)  gives 


G{4) 

<7 


sin(^T£  c)cos{0J3eL$)  -2  sin2  ^  sin (0?SEL$) 


qTE~c 


s\n{qTE+c)cos{[3jEL$)  +  2 sin2  ^  sin {0jEL%) 


qrE+c 


.  (3.108) 


Note  that  the  right-hand  side  of  (3.108)  is  that  of  (3.76)  with  <5  replaced  by 
TE  and  with  the  sign  of  the  coefficient  of  l/(qf+c)  changed.  If  we  recall  also 
that  (3.76)  reduced  to  (3.80),  we  see  that  (3.10S)  reduces  to 


+  f-l)'' 


cos(tfEL3  -  ^) 


qrr. 


(3.109) 


As  for  D in  (3.105)  and  (3.107),  substitution  of  [1,  eq.  (F.S6)]  into  (3.73) 
gives 


D™  = 


_ysin(n'r£  c)  +  2sin2 


’  TE-  ' 
n.  c 


— j  sin (nTE+c)  +  2  sin2 


'  nTE+G 


nTE~c 


+ 


nTE+c 


(3."0) 


The  right-hand  side  of  (3.110)  is  that  of  (3.81)  with  8  replaced  by  TE  and 
with  the  sign  of  the  coefficient,  of  l/(nLS-c)  changed;  recalling  t hat  (3.81) 
reduced  to  (3.82),  we  see  that  (3.110)  reduces  to 
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■  {sin  (#*£,  -  f )  +  j  cos  (f„*L3  -  f )  }  .  (3.111) 

The  F's  in  (3.105)— (3.107)  are  given  by  [1,  eqs.  (F.S4).  (F.123),  and  (F.124)] 

F (3)  =  f(nTE~c ,  - qTE~c )  +  f(nTE+c ,  qTE~c ) 

+f(nTE~c ,  <yTE+C)  +  /(nTE+c,  -r/E+c)  (3.112) 

F(1)  =  -f{nTE~c,  -qTE~c)  +  f{nTE+c,  qTE~c) 

+/(nTE-c,  qTE+c )  -  /(nTE+c,  —qTE+c)  (3.113) 

F(5)  =  f{nTE~c,  —qEE~c)  +  f(nTE+c,  r/TE-C) 

-f(nTE~c,  <yTE+c)  -  /(nTE+c,  — <yrE+c)  (3.114) 

where  /  is  given  by  (3.S4). 


3.2.4  The  TE  Quantities  for  x'rs  >  ka 

When  x'r3  >  ka,  we  proceed  to  evaluate  expressions  (3.67)-(3.70)  for  id2^  — 
—  £(Dh  Expression  (3.67)  is 


j(Fte  +  GEEDTnEhJ3Ea 


(3.115) 


where  ~JEa  is  given  by  (3.58).  In  (3.115),  j(FTE  +  GEEDEE)  is  given  by 
the  right-hand  side  of  (3.93)  with  6  replaced  by  TE,  that  is,  with  g  (which 
appears  in  the  z's  of  (3.95)-(3.9S))  given  by 


with  cn  given  by 


a  —  GTE 

y  its  i 


2~/tec 


(3.116) 


(3.117) 


n  (n~)2  -|-  (7eec)2  ’  '  ' 

and  with  c7  given  by  the  right-hand  side  of  (3.117)  with  n  replaced  by  q. 
Throughout  Section  3.2.4,  g  and  cn  are  given  by  (3.116)  and  (3.117).  This  g 
and  this  cn  are  not  to  be  confused  with  th eg  and  t.ie  cn  given  by  (3.89)  and 
(3.90)  of  Section  3.2.2. 

Expression  (3.68)  is 


i(3)  =  F(3)  +  GTqED^\ 


3.118) 
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In  (3. IIS),  F(3)  is  given  by  [1,  eqs.  (F.l  11)— (F.l  13)] 


(-l)ncnc,  sinh(jfc) 

jr( 3)  _  21  (  — I)"cncq  sinh(^c) 

5C  4{sinh(ffc)  -  gc} 

iac)2 


>  7  7^  n 

C,  ,  ?  =  7i  ^  0 

,  <7  =  n  =  0 


(3.119) 


As  given  by  (3.119),  F ^  =  0  when  q  =  n  =  0.  The  right-hand  side  of  (3.119) 
was  written  so  as  to  be  similar  to  the  right-hand  side  of  (3.SS).  In  (3.118), 
GjE  is  given  by  the  right-hand  side  of  (3.91)  with  g  given  by  (3.116)  rather 
than  by  (3.S9).  As  for  in  view  of  (3.77),  substitution  of  [1,  eqs.  (F.102) 
and  (F.103)]  into  (3.73)  gives 


n(3)  _  2n~cn  aL,  f  -  sinh(</c/2)  ,  n  even  ) 
L  "  gc  \  cosh(</c/2)  ,  n  odd  J 


(3.120) 


Comparing  the  quantities  F(3),  GEE,  and  of  this  paragraph  with  the 
quantities  Fs  of  (3.S8),  Gsq  of  (3.91),  and  Dsn  of  (3.92)  and  noting  that  the 
latter  quantities  combined  to  give  (3.93),  we  obtain 

F(3)  +  QTSfrz)  =  +  GTqEblE)}  (3.121) 

gc 

where  j(FTE  +  GEE DEE)  is  given  by  the  right-hand  side  of  (3.93)  with  S  re¬ 
placed  by  TE  as  described  in  the  third  sentence  of  Section  3.2.4.  Substitution 
of  (3.121)  into  (3.118)  gives 


Expression  (3.69)  is 

A-0  _  _(FW  +  g<4) £*•£). 

In  (3.123),  js  given  by  [1,  eqs.  ( F.  129)— (F.  131 )] 

(-l)ncnc,  sinh(#c)  ,  q  ^  n 

(  —  l)ncnc,  sinh(gc)  —  cq  ,  <7  =  n  /  0 
4{sinh(<7c)  -  gc} 


(3.122) 

(3.123) 


(3.124) 
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As  given  by  (3.12-1),  F^)  =  0  when  q  =  n  =  0.  The  right-hand  side  of  (3.124) 
was  written  so  as  to  be  similar  to  the  right-hand  side  of  (3.SS).  As  for  G^ 
in  (3.123),  substitution  of  [1,  eq.  (F.126)]  into  (3.74)  gives 

{sinh^Ig)  -  (-1)7  sinh(<7(Fj  -c))}.  (3.125) 

Substituting  (3.77)  into  (3.125)  and  using  [2,  formulas  651.06  and  651.07], 
we  obtain 


q(4)  _  :-qKcq  f  sinh(gc/2)cosh(gL3)  ,  q  even  1 
7  gc  \  s\nh(gL3)  cosh(<?c/2)  ,  q  odd  J  ‘ 

In  (3.123),  DEE  is  given  by  (3.92)  with  8  replaced  by  TE.  Thus, 

nTE  _  0  •  -gL,  J  —  sinh(<7c/2)  ,  n  even 
n  -  \  cosh(f?c/2)  f  n  odd 


(3.126) 


(3.127) 


where  g  and  cn  are  given  by  (3.116)  and  (3.117).  Comparing  the  quantities 
F of  (3.124),  6-TO  of  (3.126),  and  DEE  of  (3.127)  with  the  quantities  F 6  of 
(3.SS),  G8  of  (3.91),  and  Dsn  of  (3.92)  and  noting  that  the  latter  quantities 
combined  to  give  (3.93),  we  obtain 


fMI  +  ^  {j(Ftc  +  G^D™)} 


(3.128) 


where  j(FTE  +  GTE  DEE)  is  as  in  (3.121).  Substitution  of  (3.128)  into  (3.123) 


gives 


;(*)  =  -31{j{FTE  +  6TEt)TE)y 
gc  K  } 


Substitution  of  (3.116)  and  [1,  eq.  (F.121)]  into  (3.70)  gives 
2  a 


i(5)  _ 


1  ,  n  =  ,#0  1  (  j(F<5<  +  Gj,»P!,3>) 


cx'rl  \  0  ,  otherwise 
where  F*5*  is  given  by  [1,  eqs.  (F.  132)— ( F.  134)] 


ga 


(3.129) 


(3.130) 


F(5)  =  nfF 


(gey 


(  (-l)ncncqsinh(gc)  ,  q±n 
2  (  — l)ncnc7  sinh(gc)  ,  q  =  n  ^  0 


4{sinh(ffc)  -  gc] 

(gc)2 


+  jc, 


,  q  =  n  =  0 


1  .  q  =  n  ^  0  ] 
7  ]  0  ,  otherwise  j 

(3.131) 


22 


Note  that  [1,  eq.  ( F.  132)]  is  not  correct.  Please  correct  [1,  eq.  (F.132)]  by 
multiplying  the  denominator  of  the  right-hand  side  by  (nr)2  +  i~lJ,Ec)2-  As 
given  by  (3.131),  =  0  when  q  =  n  =  0.  The  right-hand  side  of  (3.131) 

was  written  so  as  to  be  as  much  as  possible  like  the  right-hand  side  of  (3.88). 
In  (3.130),  and  D ^  are  given  by  (3.126)  and  (3.120),  respectively.  Com¬ 
paring  the  quantities  of  (3.131),  G £0  of  (3.126),  and  of  (3.120)  with 
the  quantities  F6  of  (3.88),  Gsq  of  (3.91),  and  Dsn  of  (3.92)  and  noting  that 
the  latter  quantities  combined  to  give  (3.93),  we  obtain 


j(F<5>  +  =  2^ 


1  ,  q  =  n  ±  0  | 
0  ,  otherwise  J 


(3.132) 


where  c„  is  given  by  (3.117)  and  cnq  is  given  by  (3.94)  in  which  the  s’s  are 
given  by  (3.95)-(3.9S)  with  g  given  by  (3.116).  Substitution  of  (3.132)  into 
(3.130)  gives 


J(5)  =  m*2cnq  ,  M  (  1  >  <7  =  n^0l 

( 9a)(9c)2  Vcx'rs  9*!  l  0  >  otherwise  J 


(3.133) 
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Chapter  4 

The  Excitation  Vector 


The  elements  of  the  excitation  vector  are  given  by  [1,  eqs.  (5.15)  and  (5.16)] 

(4.1) 


where  [1,  eq.  (5.21)] 


and  [1,  eq.  (5.24)] 


JO  ,  m  even 
ysm~\  -rr  ,  m  odd 


(4.2) 


(4.3) 


sin(n7r  -  /ffi^c)  L$ )  -  2  sin2(  'c)  sin (Pff1  L$ ) 

2(n?r  -  ffi'c) 

,  sin(n7r  +  /ffiAfc)  cos (f3£{u L$ )  +  2  sin2( Lt ) 
T  2(nir  +  ${xtc) 


(4.4) 


In  (4.1)  and  (4.2),  a  is  either  1  or  2,  and  i  determines  m  and  n  as  described 
in  Appendix  A. 

Comparing  the  right-hand  side  of  (4.4)  with  that  of  (3.76).  we  see  that 


=  ;  [«"] 


oi 


(4.5) 
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where  is  G™  when  r  =  0  and  5  =  1.  Substitution  of  (4.3)  and 

(4.5)  into  (4.1)  and  (4.2)  gives 


jaTM  _  84>Qn  IT  r  ~rwi  -30™l3  f  0,  m  even  1 
'  mkmnb\l  ire  l  "  ioi  \  1,  m  odd  J 

raTE  !em^nc  f/.T.U]  .-iff™  L*  f  0?  m  even  ) 

'•  =  OVIS-^"  <01  \  1,  modd  J 


(4.6) 

(4.7) 


Multplying  both  (4.6)  and  (4.7)  by  — jeJ/3 oi'ULa  and  noting  that  em/2  =  1 
whenever  m  is  odd,  we  obtain 


-jI?™ej0™L>  = 
-jlfEe}C)™L3  = 


~j 


S4>0n 

wikmn 


•  _S 6o_  jtnc  (ATA/ 

]kmnb\2*b  \  n 


0, 

1, 

0, 

1, 


m  even 
m  odd 


m  even 
m  odd 


(4.8) 

(4.9) 
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Chapter  5 

The  Electric  Field  in  the 
Rectangular  Waveguides 


In  this  chapter,  modal  expansions  are  found  for  the  electric  fields  in  the 
rectangular  waveguides.  Afterwards,  the  time-average  powers  of  the  TE\q 
modes  in  the  rectangular  waveguides  are  obtained. 


5.1  Expansions  in  Terms  of  the  Fields  of  the 
Magnetic  Currents 

The  electric  field  in  the  left-hand  rectangular  waveguide  (region  1)  in 
Fig.  1.2  is  due  to  M ^  and  is  given  by  [l,  eqs.  (2.3)  and  (2.11)] 


£"'  =  £  E  C"£("(!i,m;™)  +  E 

7= 1 p= 1  7=0 


E  CE£"'(a,ii^£).  (5.i) 

p=0 

p+q^tO 


In  obtaining  (5.1),  the  upper  limits  on  p  and  q  in  [1,  eq.  (2.11)]  were  sup¬ 
pressed.  We  truncated  the  double  summations  in  (5.1)  by  retaining  only 
terms  for  which  both  p  and  q  are  so  small  that,  according  to  (A. 2), 


(Ptt)2  +  <  BKM. 


(5.2) 
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The  electric  field  E}2^  in  the  right-hand  rectangular  waveguide  (region  2)  in 
Fig.  1.2]  is  similarly  given  by 

£,!|  =  EEr'£'”(a +  Z  Z  (5.3) 

q=  1  p=l  q=Q  p — 0 

P+7#0 

We  approximate  the  M’s  on  the  right-hand  sides  of  (5.1)  and  (5.3)  bv  the 
M’s  given  by  [1,  e*qs.  (3.3)  and  3.13)].  In  Sections  5.2  to  5.5,  we  express  the 
resulting  approximate  E’s  on  the  right-hand  sides  of  (5.1)  and  (5.3)  in  terms 
of  the  modes  of  the  rectangular  waveguides. 


5.2  The  Electric  Field  of  the  Magnetic  Cur- 

-  1 TM 

rent  M, 


-pq 


The  electric  field  £(1)(Q,  XI  pq  )  due  to  XI  pq  in  region  1  of  Fig.  1.2  is  given 

by  [1,  eq.  (3.30)] 


£"'(0.A/™)  =  k”V,--+)  - 

Ipq 

From  [1,  eqs.  (A. 3)  and  (A.  13)],  the  mode  field  EPgf~  is  given  by 
E™~  =  -Z™{£”(y+,z+)  - 

Ipq 

In  view  of  (5.5),  we  recast  (5.4)  as 


(5.4) 


(5.5) 


(5.6) 


5.3  The  Electric  Field  of  the  Magnetic  Cur- 
-2  TM 
rent  Mm 

m  2TM  2TM 

The  electric  field  £^(0.  Mpq  )  due  to  M  pq  in  region  2  of  Fig.  1.2  is  given 
by  [1,  eq.  (3.36)] 

Zm(b,iC")  =  (4?V,t+)  +  (5.7) 

(pq 
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(5.8) 


From  [1,  eqs.  (A. 2)  and  (A.3)j,  the  mode  field  Kpgf+  is  given  by 
E™+  =  Z™{tZU(y+,z+)  ~—}e~ 


Ipq 


In  view  of  (5.S),  we  recast  (5.7)  as 

(5.9) 

5.4  The  Electric  field  of  the  Magnetic  Cur- 

~1TE 

rent  Mp, 


The  electric  field  £(1*(0.  A/pj  )  due  to  A/p?  in  region  1  of  Fig.  1.2  is  given 
by  [1,  eqs.  (3.32)  and  (3.34)] 


)  =  iEV\*+) 

jain{ML,  V  J  1 1  +  ZiF^Eco$(/j1o(Z.i  +  r)) 
j  sin(^i0a:1)  +  ZlY1lE  cos(/^i0.ri) 

£(1)(2„&0  =  £E(si+,»+)e'*'l’+"\  (p,<l)?(  1,0)  (5.11) 

where 

Xi  —  Li  ~  xQ.  (5.12) 

In  obtaining  (5.10),  we  substituted  jf3 10  for  710  and  used  [2,  formulas  654.6 
and  654.7],  Equation  (5.10)  is  recast  as 


£("(Q,m!oE)  = 


(ZXY™  +  l)eJfr<><L«+*>  4-  {ZXY&E  -  l)e-^lo(L|+r) 

2 {j  sin(/3i0xj)  +  Z 1  Yeqe  cos( ,di0.r  1 ) ) 

•ef(f('/+>2+)-  (5-13) 


From  [1,  eqs.  (A. 14)  and  f A.  15)] ,  the  mode  fields  Ejo+  and  Ej0E  are  given 

by 

£[oE+  =  eJoE(y+,z+)e-^°x  (5.14) 

Elo~  =  fifffoVV*0*,  (p.q)  ?  (1,0).  (5.15) 
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In  obtaining  (5.14)  and  (5.15),  we  substituted  j3i  for  710.  In  view  of  (5.14) 
and  (5.15),  we  recast  (5.13)  as 

£«‘)(0  M'TE\  =  (Z-r.oE  +  eja,Dt|  +  (ZiKf  - 

10  2(j  sin(/?i0Xi )  +  ZiYwB  cos(fliox,)) 

(5.16) 

From  [1,  eq.  (A.  15)],  the  mode  field  E_ ^  is  given  by 

RTpqE~  =4qE(y\z+)e^.  (5.17) 

In  view  of  (5.17),  we  recast  (5.11)  as 

(p,?)/(l,0).  (5.18) 

5.5  The  Electric  Field  of  the  Magnetic  Cur- 

-2  TE 

rent  M;l, 


The  electric  field  £<2)(Q,iL Lpq  )  due  to  \Lpq  in  region  2  of  Fig.  1.2  is  given 
by  [1,  eqs.  (3.38)  and  (3.40)] 

r(2)m  m2TE\  -  J  sin(/3io(L2  -  x))  +  Z2Y^E  cos{(3l0{L2  -  x))  TE  +  + 

£l  (Li,  M-io  )  \  T  7  vsTE  (a  \  —10  ( V  iz  ) 

j  s\n(f3mx2)  +  Z2  Vio  cos(/?,02-2) 

(5.19) 

£(2,(Q„0  =  £E(!/V-->-W'-">,  (m)#(1,0)  (5.20) 

where 

x2  =  L2-  Xo.  (5.21) 

In  obtaining  (5.19),  we  substituted  j/3l0  for  710  and  used  [2.  formulas  654.6 
and  654.7].  Equation  (5.19)  is  recast  as 


£™(Q,M;f) 


(ZiY™  +  +  (Z2Y]lE  -  l)e-2^o(^-J) 

2 (j  sin (/3lox2)  +  Z2\\ T0E  cos(310x2 )) 

(5.22) 
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In  view  of  (5.14)  and  (5.15),  we  recast  (5.22)  as 

pmm  {Z2Y™  +  1)E™+C^>  +  (Z£gf  -  l)£j0E-e-is'°L> 

10  2(j  sin(/3ioX2)  +  22>h(£cos(/?i0i2)) 


From  [1,  eq.  (A.  14)],  the  mode  field  £p,E+  is  given  by 

£jf+  =  <'V,z+)e"’"1. 

In  view  of  (5.24),  we  recast  (5.20)  as 

»  2T  E 


Em(a,Utt  )=ElE+e^,  (p,q)?(  1,0). 


(5.23) 


(5.24) 


(5.25) 


5.6  Expansions  in  Terms  of  Waveguide  Modes 

Substitution  of  (5.6),  (5.16)  and  (5.18)  into  (5.1)  with  the  AT s  replaced  by 
A/’s  gives 

,-ji3\oLi 


Em  =  (Z,K,7  +  +  (ZiYtf  ~  ')£??+ e 

—  21  j  sin(/)10n)  +  Z\ >',oE  cos(Aon)) 

7=1  P=1 

+  E  E  «£'eWt. 

q=0  p=0 

p+7/0 

(pw)#(i.o) 


i/lTE 
"  V  10 


(5.26) 


Substitution  of  (5.9),  (5.23)  and  (5.25)  into  (5.3)  with  the  A/’s  replaced  by 
A/’s  gives 

£m  =  ± i)£foE^,l,10tl  +iM° E  -  i)£f0E^~rf’ofa  r«-E 


2{j  sin(/?ioi2)  +  Z2V12e  cos(/?,0x2}) 

+LEC',(4C)'w' 

P'7 


10 


7=1  p=l 


+  E  E  v™  %?*<?’'”■ 


(5.27) 


q=0  p=0 

p+l/o 

(p.q)^(l,0) 
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5.6.1  Normalization  of  £(1^  and  £(2) 

The  quantities  E and  E are  due  to  the  2  traveling  wave  whose  electric 
field  in  the  circular  waveguide  is  ikJiUe+  given  by 


=  zir°eir^9)i 


(lcTAr 12 1 bT A/e ( o  6)e -i*3™ 2 

“o»  -  +Mzizoi— Liioi  -  (5.28) 


where 


yTAfeo 

^01 


'/A 


TA/ 

01 


(5.29) 


Equation  (5.2S)  was  obtained  by  substituting  (3.57)  into  [1,  eq.  (B.l)].  Equa¬ 
tion  (5.29)  was  obtained  by  substituting  (3.57)  into  [1,  eq.  ( B .25)] .  In  this 
subsection,  suitable  expressions  are  found  for  the  quantities  E(1  )ejPoistL3  / 

\J Z0™«  and  E}2'1  eJ Z^,eo.  These  quantities  are,  according  to  (5.28), 

due  to  the  2-traveling  wave  whose  transverse  electric  field  is  \J Zj'1UeoeJA/  e(/>» 
at  2  =  L3  in  the  circular  waveguide.  In  Section  6.3,  it  will  be  shown  that  the 
2-directed  time-average  power  associated  with  this  field  is  unity. 

Multiplying  both  sides  of  (5.26)  and  (5.27)  by  eJl3°*'L3 / \J Z2)y!ro  and  using 
(5.29),  we  obtain 


Mi 


Meo 


r<\TE- piPxoio  pTE-  ,  p-lTE+  -j0iOx„  pTE+ 
1)  _  °io  e _ ^\o_  ^  L  io  e  £-10 


TE 

10 


r\TM- eipqx0pTM- 

V'  V'  — El  — P7  ST 

/l  ?T  \/  1  ' 

q=ip=i  y'^pq  I  7=0  p=° 

p+7^0 

(p,q)#(l,0) 


r\TE-  „7p,x„  pTE- 

C'p7  c  ^-pq 


# 


T  E I 
pq  1 


(5.30) 


.y/Z{ 


T  Meo 
01 


£(2)  = 


C?IE+ eJl3'aXo  ET,?+  +  C2[E-e-^lor°£[0E 


V7 


-TE 

10 


p2TM+  -tpqx0  pTM+ 

+  EE  ”■  >  +  E  E 


q=l p=l 


7=0 


p=0 

P+7^0 

(p,7)^(1.0) 


p2T E+  p~)PqX0  pTE+ 

- _^7— .  (5.31) 

■'TCI  ' 
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V-2TEeJJ0V'Z.3 


(5.37) 


cm,+  _  jk 

’ "  r  \  i  ) 


r‘2TE+ 
^ P9 


y2TEe:3^L3  \ 

*  J 


(5.38) 


(5.39) 


5.7  Time- Average  Power 


What  is  the  x-directed  time-average  power  associated  with  an  arbitrary  elec¬ 
tromagnetic  field  (£,/?)  in  a  source-free  region  of  a  rectangular  waveguide? 
This  field  can  be  expressed  as 


E  =  E(C™+£”'+  +  +  Cjf-i™- )  (5.40) 

P.9 

a.  = 

P.9 


iTM+  ltTM+ 


iT\f-  uTM- 


-  TTl  M+  ,  ^1  M  —  rrl 
'  P9  —P9  ‘r  —P9 


+  c™+a™*  +  c™-h£- ).  (5.4i) 


-P9 


P9 


-P  9 


On  the  right-hand  sides  of  (5.40)  and  (5.41),  the  C’s  are  constants,  and  the 
Ks  and  the  /£’ s  are  the  mode  fields  defined  by  [1,  eqs.  (A. 2).  (A. 3),  (A. 14), 
and  (A. 15)].  The  x-directed  time-average  power  P  associated  with  (£,//)  is 
given  by  [4,  eqs.  (1-57)  and  (1-58)] 

P  =  [  dy+  [C  dz+Re(E  x  H‘)  •  ux  (5.42) 

Jo  Jo 

where  denotes  complex  conjugate  and  uRe'’  denotes  real  part.  Sub¬ 
stituting  (5.40)  and  (5.41)  into  (5.42)  and  using  the  previously  mentioned 
definitions  of  the  mode  fields  and  the  last  of  the  orthogonality  relations  [1, 
eq.  (A. 26)],  we  obtain 

P  =  -  \C™-\2e2n'^}Rc(Zj") 

P.9 

-2  E  Imag{C™+(C™-  )-e-^'"“*<’«>')Imag(  Z™ ) 

P.9 
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In  (5.30) 


r\TE- 

Wo 


(Z.Y^  +  ^e 


\ 


2(ZiVjo£  cos(/310-Ti )  +  j  sin(/3i0xi))  J 


(5.32) 


C}TE+ 


10 


I  P 10 

fiTU 

Poi  , 


(Ziffi-  l)e-^10J- 


2(ZiVjjE  co*(3 io.i*i )  +  j  sin(/?ioJ’i ))  J 


y/lTE  jJqi  Lj  \ 

1 - )  (5.33) 


C, 


\TM-  _ 
Pi 


c. 


ITE-  _ 

P9 


jfe 

y\TE Ls  \ 

DO  I 


n 


4Af 


(5.34) 


(5.35) 


In  (5.31), 


r>2T£+ 

Wo 


I  /3l0 

qTM 
Poi  , 


(Z2y: 


-T£+;)eWoo 


2(Z2V^£cos(^o.r2)  +  7  sin( 3l0x2)) 


i/2 TE  }3£"U\ 

ii2-l -  (5.36) 


C 


2TE+ 

10 


/  HW~\  f  (Z2^c-  De-J^ 

\V  /3m U  j  \2(Z2V'17u£'co^  iio.1'2)  +7  sin(  J,o.r2)) 
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+  E{!^E+|2e"2Re<"M)X  -  |CpT/-|2e2Re,-'-)'r}Re(Vp7) 

p.? 

-2^Imag{CXEt(C™-)-e-2J|"'**<^''Imag(V^£)  (5.43) 


where  “Imag"1  denotes  imaginary  part.  In  (5.43),  we  have  [1,  eqs.  (A. 13)  and 
(A. 25)] 


7TM  _  :  VV pq 

"  " 

yTE  _  ■  7 pq 

pq  J  ku 


(5.44) 

(5.45) 


where  7 pq  =  —jf3pq  if  the  mode  propagates,  and  ~fpq  is  purely  leal  if  the  mode 
does  not  propagate.  Since  only  the  T Ew  mode  propagates,  (5.43)  reduces  to 

P  =  J  {EW.nag{C™*(C” •'-)■}}  +  ^<ICroE+l2  -  \CLE-f) 


+ir  E7„'mag{cEE+(cPT,E')')  • 


(5.46) 


c<? 

(p, 7)^(1, 0) 


5.7.1  Time- Average  Power  in  the  Rectangular  Wave¬ 
guides 

The  normalized  electric  fields  £}  1'eJ/3oi'U L* / yj Zjxsteo  and  E}2^e3^™ 1,3  / 
\J of  (5.37)  and  (5.3S)  are  due  to  the  ^-traveling  wave  whose  trans¬ 
verse  electric  field  is  E^^re+c"3f}^'  ,^Z~L^/\J Zq j-Ueo  in  the  circular  waveguide. 
The  r-directed  time-average  power  of  this  field  in  the  circular  waveguide  is, 
as  given  by  an  expression  very  similar  to  (5.31),  equal  to  unity.  The  — x- 

directed  time-average  uower  of  gj^e313™  Ls  / \J Zjj'Ueo  is  P^/Z™1*0  given  by 
(5.31 )  as 

=  |Ci'oTE'l2  -  |c;„T£+p.  (5.47) 


34 


The  ;r-directed  time-average  power  of  / yj Z™eo  is  pP  /  Z™eo 

given  by  (5.31)  as 

=  |C?0re+|2  -  |C?„re-|2.  (5.48) 

^01 

When  the  incident  time-average  power  in  the  circular  waveguide  is  unity,  the 
time-average  power  Pt  transmitted  into  the  rectangular  waveguides  is  the 
sum  of  (5.47)  and  (5.48): 

Pt  =  ICS*- 1*  -  ICTE+|2  +  !C,2oT£+l2  -  |C?0TE-|2.  (5.49) 
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Chapter  6 


The  Electric  Field  in  the 
Circular  Waveguide 


In  this  chapter,  a  modal  expansion  is  found  for  the  electric  field  £(3*  in  the 
part  of  the  circular  waveguide  for  which  zs  <  z  <  —c/2  where  is  such  that 
the  impressed  source  1Z.imp  lies  in  the  region  for  which  z  <  zs.  See  Fig.  1.2 
where  the  entire  region  inside  the  circular  waveguide  is  called  Region  3.  As 
stated  in  Chapter  1,  only  the  TM0 1  and  TEn  modes  propagate  in  Region 
3.  The  coefficients  of  the  TM0 j  and  T En  modes  in  the  modal  expansion 
for  are  then  expressed  in  forms  suitable  for  computation.  Finally,  the 
time-average  powers  of  the  T Mq\  and  TE\\  modal  contributions  to  are 
obtained. 

The  electric  field  in  the  circular  waveguide  is  given  by  [1,  eq.  (2.7)] 

E(3)  =  E(3)(0,  — A/(l)  -  A/(2))  +  £(3)Uimp,Q)  (6.1) 

where  E(3)(Jimp,  0)  is  the  electric  field  due  to  Jimp,  and  E(3)(0.  -3/(1) -il(2)) 
is  the  electric  field  due  to  —M}1'1  —M}2^  where  A//1'  —  is  the  combination 
of  — and  — A/(2).  Each  of  the  sources  JZ’mp,  —  ,U(1),  and  — radiates 
in  the  circular  waveguide  with  the  apertures  closed,  with  the  short  at  z  =  L3, 
and  with  a  matched  load  at  the  other  end  where  z  <<  0. 
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6.1  The  Electric  Field  £<3)(Q,  -M(1)  -  M(2)) 

Since  M(1)  and  are  given  by  [1,  eqs.  (2.11)  and  (2.12)],  we  have 

£(3)(Q,-M(1>  -m(2))  =  -  £  V;™£(3)(Q,i/2™) 

-  £  V7,te£(3)(Q,M^te).  (6.2) 

The  first  12~,,q,P  on  the  right-hand  side  of  (6.2)  stands  for  Z^=i  H,=i  ZP=i- 
The  second  Yl-y,q,P  stands  for  Ylq=o  Z  p= o  •  Here,  no  upper  limits  are 

P+<7#0 

placed  on  the  indices  p  and  q.  Terms  are  retained  only  for  which  p  and  q 
are  so  small  that  (5.2)  is  true.  The  ’s  on  the  right-hand  side  of  (6.2)  are 
given  by  [1,  eqs.  (4.47)  and  (4.61)]: 

Fmln  .,,™> _ “a  .  q  ~  js.  f.(e,)v:(g,°)gT-w 


(kJffrJr(k™aWTEt 
c  ke«)(k2,  -  r2)Jr(K,) 

,  p ' ,  sin  A, ,  j?.  e, ( 4f,E)3  Jr  ( ■ kjfa )-FTE-'  1 

(t,(  A,  ’.tsfe  2-,™(I;5  -  r»)^(4,)  f  (  '  > 


£,3’(a  .*£*)  = 


4a  ,p  “  “  e,(*™)V;(*™0)£>™* 


+(f)ES 


4  .tss;  2x;,j?+1(ir,) 

(klEfrUk™a)FTBi 


+(c)(— (6-4) 


In  (6.3)  and  (6.4),  is  given  by  (3.58).  Furthermore,  £7rA/i',  E?TE'b ,  and 
£PTEz  are  given  by  [1,  eqs.  (4.74)-(4.76)]: 


7**7  T  M<t> 


=  {-;™1  sinh(7rT;w(L3  -  *))  + 

f  Vllr™  n\  -  07lr-4(^'VV)  I 

•  <  Mp<?  JT\^rs  P )  ^  kT^  p  J 


^v:(fcr?v)-i4- 
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+  {*™  cosh(7™(£3  -  .-))  +  af,k™  if™'*  (6.5) 

5Vi 

t?TB*  =  {-2TEIsinh(7;£(£3-J))  +  ireV’”|1>-'} 

'  -  ***"■>'(&*?)}  <6-«) 
JjpTEz  __  ^ZTE3  s[n\i^JsE^ll3  ^  z))  +  ZTE‘ie~'y^E('L3~z^ 

■(-i  r  -  }  ■  (W 

The  <f>' s  in  (6.5)-(6.7)  are  given  by  [1,  eqs.  (E.10)-(E.13)]: 


=  -(-ir{ 


41)cos(— — )  -  42)  sin(  — — ) 


4?2  =  <f>(2)  cos(  — — )  +  4>W  sin(— — ) 

X  Q  X  Q 

4>l3  =  -(-IT  |<^3)cos(  — — )  -  <^i4*  sin(  — — ) 


(6.10) 


4>lA  =  ^  cos(— — )  +  43)  sin(— — ). 


(6.11) 


The  y's  on  the  right-hand  sides  of  (6.S)-(6.11)  are  given  by  [1,  eqs.  (2.15) 
and  (2.16)]: 


y1+  =  (ir  -  <t)z,  +  ~ 


(6.12) 


2+  i  u 

y  T  =  <px0  + 


(6.13) 


The  <f>p' s  on  the  right-hand  sides  of  (6.8)-(6.11)  are  given  by  (3.40)-(3.43). 
For  simplicity,  we  assume  that  z  <  —c/2  so  that  [1,  eqs.  (4.77)  and  (4.83)] 


ZTM\  _  ZTE\  __  JTEZ  _  Q 


(6.14) 


Since  z  <  —c/2,  the  remaining  superscripted  z’s  in  (6.5)-(6.7)  (namely  zr‘U2, 
zTE\  and  zTE4)  are,  as  stated  in  (1,  page  34],  given  by  their  expressions  in 
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[1,  Appendix  F]  with  z+  replaced  by  zero.  For  xrs  <  ka,  we  have  [1,  eqs. 
(F.26)  and  (F.35)] 


f!! 

c 

,  + 

jtea 


|sii 


sin(<75  c)cos(£sL$) -2sin2(^ -)sm(f3srsL$) }>  /(2q5  c) 


2)} 


!  sin(/+c)  cos (f3sraL$ )  +  2  sin2(~)  sin(/3^Z^)|  /(2 qs+c) 


=  -j 


(6.15) 


f  TE —  A 

sin (qTE~c)  cos (/3jsEL+)  -  2s\n2{~-^)sm(/3jE  L^)\  l{2qTE~c) 


+j 


sin {qTE+c)  cos(/3rr/l+ )  +  2  sin2(  sin(/3rT/L+)|  /(29T£+c) 


(6.16) 


For  xT,  >  lea,  we  have  [1,  eqs.  (F.33)  and  (F.42)] 

2 
A. 

c 


.TEA 


=  {^.csinh^ij )  -  —  c))}  /  {l'/«  )2  +  (l,*,c)2} 

(6.17) 

=  <jir  {sinli(7r7,El.J)  -  (-1)’  sinh(7-T,E(/.J  -  c))}  /  { f .J-) ’  +  (7rTEc)2} . 

(6.18) 


In  (6.15)  and  (6.17),  S  is  either  TM  or  TE. 

'  6.1.1  The  Quantities  £rW,  EyTE4> ,  and  £7TE‘~ 

In  this  section,  expressions  (6.5)-(6.7)  for  _E'7TA/<®,  jpTE* ,  an(j  jjnTEz  are 
reduced  by  means  of  (6.14),  and  then  expanded  by  means  of  (6.S)-(6.13). 
Next,  the  —2-traveling  modes  of  the  circular  waveguide  are  introduced.  Fi¬ 
nally,  the  expressions  for  EP™^ ,  EPTE<t’ ,  an(j  ePTEz  that  were  obtained  by 
means  of  (6.8)— (6. 14 )  are  recast  in  terms  of  the  -2- traveling  modes  of  the 
circular  waveguide. 
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Reduction  of  (6.5)-(6.7)  by  means  of  (6.14) 

Substitution  of  (6.14)  into  (6.5)-(6.7)  gives 


f?~<T M<t>  _  zTM2e-y™{L3~z) 


f.,  --V2  7//.TA/  X  ,  <PlrJT(k™p)  4>y2k™Jr{kJsMp) 

Jr(kT3  p)-  24 - jjjj- - +  MZ - - 


z™  = 

£',TEi  =  ^£'>e-’ie(t«->)(_i)-r  L 

l  P 


(6.19) 


(6.20) 


+  U<j><t>'l4J'r[kJsEp)  >  .(6.21) 


Expansion  of  (6.19)-(6.21)  by  means  of  (6.8)-(6.13) 

Substituting  (6.12)  and  (6.13)  into  the  arguments  of  the  trigonometric  func¬ 
tions  in  (6.8)-(6.11)  and  using  [2,  formulas  401.01-401.04],  we  obtain 

C°S  =  (-irr  (cos  cos(r^)  -  (  — l)‘v  sin  sin(rd>)  1(6.22) 


=  (  — l)’1"’  s  sin 


cos(r<£)  -|-  (  —  l)7  cos 


Substitution  of  (6.22)  and  (6.23)  into  (6.8)-(6.11)  gives 


sin(r^)  >(6.23) 


f1  =  (~l)7r  {(-1)7+VJ1  cos(r<?i>)  +  <t>bp2  sin(r^)} 
<P2  =  (-lTr  cos(r^)  +  (-l)^J1  sin(r^)} 

=  (-irr  {(-l)7+1^3cos(rd>)  +  sin(r^)} 
<P4  =  (-ir{^4cos(r^)  +  (-l)^3sin(r^)} 


where 


^=^1)C0S  (^)'^Sin(^) 

€  =  0<2)cos  +  ^’sin 


(6.24) 

(6.25) 

(6.26) 

(6.27) 

(6.28) 
(6.29) 
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4>f  =  4 >i3)  cos  ( -  4>l4)  si 


2x0 


sin 


rb 

2Z 


=  ^’cos  (^)  +  43,sm  (^-j 
Substituting  (6.24)-(6.27)  into  (6.19)— (6.21),  we  obtain 

jjpTMtt>  _  ^\yr2™2e-‘ $ 


(6.30) 

(6.31) 


£ 


sinW)+!4«^M 

k™JT{k™p)  sm(r4>) ,  ji>2f  t/ / »  t\[  \  ,  n 

+uz - ~jn{7 - }  +  <Pp  {MpJr(tr,  P )  cos(r^) 

•  rs 

}}  (6.32) 


-14 


rJr{k™p)sin(r4>)  k™  JT{k™ p)  cos(r<£) 
- (-  U, - =rr7 - 


£-yTE<t>  _  2TE2e—yJ'sE{L3-z) 

(-irtfK^(^g;-B(^  +  ~W)}  - 

+^Mng—  - 14 J'AkJfp) sin(r^)} }  (6.33) 

jriTEz  _  {—\yrzTE*e-' (£*-*)  ■£ 

-tfKrJr^”(^  -l^erisinW)}  +  (_ir 

,64  f  rjr{kjfp)sm{r4>)  Te  s  ,  \ 

■4>„  {lip - - +14Jr(fcr,  /?)  cos(r4>)  j  j  .  (6.34) 

The  -2-Traveling  Modes  of  the  Circular  Waveguide 

The  modes  of  the  circular  waveguide  that  travel  in  the  -2-direction  are 
£™",  £lU°-,  £,E,~,  and  ZlE°-  given  by  [1,  eqs.  (B.2),  (B.27),  (B.36), 
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and  (B.56)] 


=  | ~Z™eoe™e(p,  <f>)  +  uz |  e^M*  (6.35) 


?T\fo- 


_  )  yT.\leo„T\lor  i\  .  ,, 

=  )  ~^rs  &r,  (Pi  <P)  +  lLz 


(6.36) 


KlEe~  =  £*(p,*)e**  (6.37) 

Ej?°~  =eJsEo(p,6)e*’E*  (6.38) 

where  [1,  eqs.  (B.25),  (B.7),  (B.22),  (B.30),  (B.33),  (B.ol),  and  (B.62)j 


JTM 

7TXIeo  _  its 


V  7T  3?ra «7r+l(^ra) 


£r?">,<?)  = 


wrM)  =  \  - 


7T  \ajr+1(xr3) 


•  cos(ro)  -  it, 

[2Jn(k^p)sm(n<i>) 


rJT(k™  p)sm(rp) 
k™P 


ft  ■Enp'Ai+l  (^np) 


£A"W)  =  -1/3 


ft  \Q>/r+l(^'ri) 


U-pJrik™ p)  sin(r<£)  +  it, 


rJr(kJsAIp)  cos(rd) 

’  ev 


£*•(/».  *)  = 


7r(x'r2-r2)  \  ./,(*'„) 
rJr{kEE p)  sin(rei>) 
"p  *£*/» 


+  ILi^(^r/p)cos(r0)  1  (6.44) 


(6.39) 


(6.40) 


(6.41) 


(6.42) 


(6.43) 


ft(x?s  ~  r2)  \Jr(Xrs) 
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.  +  lL^J'T{kJEp)  sin(r^)|  .  (6.45) 

In  view  of  (6.39)  and  (3.49),  substitution  of  (6.40)  and  (6.41)  into  (6.35)  and 
(6.36)  gives 


?TMe-  _  fir  {  Zj,X! 

V  7T  ^  ^ 

rJr(k™p)  sm(r4>)  |  ^  k™JT{k™p)  cos(r^)  1  ^  ^ 


?TMo- 


<±4 - 

2  /  Z™ 


Ar™p 


:  {l±pJ'r(k™ p)  cos(r<f>) 

+  1±, 


rs 

to 


<yTM 

ITS 


rs  _  1  e7„* 


+lhf 1 


,\^+,(x„);-™={M^'”V)sinW) 

rjr(k™p)  cos(rd)  k™Jr{k™p)  sin(rd)  \ 


+  2i* 


'  “*  7™  J 

Substituting  (6.44)  and  (6.45)  into  (6.37)  and  (6.38),  we  obtain 


.  (6.47) 


?TEe— 


e**. 


“  r2) 

(6-48) 


?TEo- 


^  \  e~iT*z 


n(x'r]-r2)  \Jr{x'r ,) 
rJr(kffp)  cos(r<f>)  TE 


lip 


kJsEP 


+  lLfj>J'r{kr3  p)  sin (r<£)  >  .  (6.49) 


Expressions  for  £7™^,  EyTE4>,  and  ElTEz  in  Terms  of  Waveguide 
Modes 

Equations  (6.46)  and  (6.47)  reduce  (6.32)  to 


z 


TM2 


t.  \  z™’° )  aJ't,U”'1 


(6.50) 
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Equations  (6.4S)  and  (6.49)  reduce  (6.33)  and  (6.34)  to 

£>™  =  (-1  -  r2) 

•  {(-ir<#£^e-  _  ^£T£o-  j  e-,TEL3 

r™  -  (-i 


{(-ir^,£t-  +  <3£"“-}e 


—yJ.EL3 


(6.51) 


(6.52) 


6.1.2  Expression  for  £^3)(Q,  —  —  A/^)  in  Terms  of 

Waveguide  Modes 

In  this  section,  expressions  (6.50)-(6.52)  for  £'l'TiU'*>.  E^TE6 ,  and  EPTEz  are 
substituted  into  expressions  (6.3)  and  (6.4)  for  E(3)(0.  Mp™1)  and  £(3)(0, 
M£E).  The  resulting  expressions  for  £(3)(Q,  ALpJ*1)  and  E{3){0.,  APpJE)  are 
then  substituted  into  expression  (6.2)  for  £(3)(Q,  — J/(1)  —  AI(2)). 

Substituting  (6.50)-(6.52)  into  (6.3)  and  (6.4)  and  using  (3.49)  (3.51), 
we  obtain 


cZTMe° 


9  lOirh  f  00  /T —  00  XA/2 

£™(Q,A/r)  = 

KP1°  '  L  r=0  “  3=1  CZ/>3 

oo  oo  yTE2 

■  «J£™'  +  (-1)"<£”'°-)  -  9 E(-*rr  Z  7=  - 

r=l  s=l  C\Jxrs  r 

•  -  tftZ?’-)  e-1”11  +  f  (^)  D-ir  yf 

■  £)  — ((-l)V“£i&-  +  o”£r,£”-)  (6.53) 


oo  „X  E2 


3=1  cJx'r2  -  r2 


J2JTE\ 


3=1  ~iJEacJx'r 2  -  r2 


2  2t rb  [e^fpc 


oo  .T.U2 


£,3,(“’^T£)  -  i^VcV  4 


•  (<A“£™'  +  (-ir<‘£™“-)  -  'j  E(-l)’r.-  £ 


oo  _T  £2 


3=1  C\  x'r\  -  r2 
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In  (6.61)  and  (6.62),  =  0  when  p  =  0  or  q  =  0. 

6.2  The  Electric  Field  £(3)Uimp,Q) 

Since  ar01  <  ka,  we  have,  from  (3.57), 

7i Tn"  =3^.  (6-63) 

Substitution  of  (6.63)  into  [1,  eq.  (5.10)]  gives 

rtry)  =  -4 — r7-T{ttA™^.(CV)sin(A™(i,-.-)) 

\yKQU!€.Ji[XQi  J 

-lL,k™J0(k™p) cos -z))}  (6.64) 

where  (3™  is  given  by  (3.59). 

6.2.1  Expression  for  £^)(  J""p,  0)  in  Terms  of 
Waveguide  Modes 


Since,  as  stated  in  [1,  page  37],  the  only  r-traveling  wave  contained  in 
Ri3)U'mp,Q)  is  the  unit  amplitude  z  traveling  T A/0i  wave,  we  suspect  that 
£(3)(i'mp ,Q)  contains  only  T M0i  modes.  Expression  (6.64)  for  ,Q) 

is  recast  as 


fifUt t4p) 


+  1L . 


we 


JU-’C 


}«- 


rM , 
i  * 
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,-2j^3  r  ,  „  k^Mk^p)y^tUi 


{  -Jtp 


+  ILz 


\ZraJi(xox)  v  "  we  ~  jwe 

From  (6.40)  and  (6.41),  we  have,  in  view  of  (3.46), 

■WffV) 

\/7rxoiJi(xoi) 

Wp) 


eJl3°' 


(6.65) 


4>)  = 


£?"'(/>»  =  -11, 


y/ffaJi  (xoi) 

so  that  (6.65)  becomes,  with  the  help  of  (3.51)  and  (3.49), 


(6.66) 

(6.67) 


}«-■ 


{fli  JVl 

we  jue 

+e-w.',"i »{  _  +  (*y)Vg»v»> },«;»■«.  (6 

^  we  jcwe  J 


.68) 


Substitution  of  (6.63)  into  (6.39)  gives 

A™ 


01  _  yTMeo 


we 


'01 


Substituting  (6.69)  into  (6.6S)  and  recalling  (6.63),  we  obtain 

£(3)Uimp,Q)  =  £™e+  4-  e-2j#*M£™e- 

where  Koife+  and  £j(Ue_  are  given  by  [1,  eqs.  (B.l)  and  (B.2)]. 

6.3  The  Electric  Field 

Since  (6.47)  gives 

£jivr°~  =  0 

and  since  (3.57)  gives  (6.63)  and 

In  =  j/Jff5, 


(6.69) 


(6.70) 


(6.71) 

(6.72) 
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substitution  of  (6.55)  and  (6.70)  into  (6.1)  yields 


(3)  _  pTMe+ 


E ™  =  Eii 


+  2 


^ bS%( 


TMe.-]0™L3 


ZTMeo 
01 


+  |  £™« 


+ +  SS0^)  +  2\j-~  £  ^ 

{  E 


5=1 

(r,j)^(0,l) 


ZTUeo 


+  E  tf,E'tZ,E’-  +  SlE°E™°-)e- 

5=1 


yJ/E: 


} 


(6.73) 


The  only  --traveling  wave  on  the  right-hand  side  of  (6.73)  is  the  incident 
wave  whose  electric  field  Ej,Ue+  is  given  by  [1,  eq.  (B.l)] 

pTMe+  __  f  rjTMeo  T  Me  /  i\  ,  (^u'OVoi'N/h  &)  ]  3 T"  z 

£01  =  )  An  £01  (P,<?)  +  lLz - : - fe  01  (6.i4) 

The  time-average  ^-directed  power  associated  with  Eo\,e+  is.  as  given  by 
a  formula  similar  to  (5.31),  Z™*0.  Moreover,  £j1'/e+  has  the  phase  factor 
e-j0oistL3  when  c  =  L3.  However,  the  normalized  electric  field  E^\,e+  eJ‘3™ Li  / 

\fzZE°  has  unit  power  and  no  phase  factor  when  ~  =  Z,3.  Thus,  the  z- 

traveling  wave  part  of  /yf Z™'0  has  unit  power  and  no  phase 

factor  when  z  =  L3.  Multiplying  both  sides  of  (6.73)  by  e*13™ 13 / \J Zj*Ieo, 

we  obtain  the  normalized  electric  field  Li  I  \jZQ\Ieo  given  by 

/  iM^3 


sJW' 
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e30ilML3 
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oo 


+  £ 


£ 

3  =  1 

(r, «)#(!,!) 


+  c„ 


TEo  rpTEo— 
rs  “rs 


) 


(6  75) 


where  Z™eo  is  given  by  (6.39)  and  [1,  eq.  (B.54)] 


yrTEeo 
*  rs 


JTE 

Irs 


JUJfi 


(6.76) 


In  (6.39),  7™  is  to  be  replaced  by  j/3™.  In  (6.76),  is  to  be  replaced 
by  jflxf.  The  C’s  in  (6. 7  j)  are  given  by 


(6.77) 

(6.78) 

(6.79) 

(6.80) 

(6.81) 

(6.82) 

(6.83) 


The  magnitudes  of  the  squares  of  the  constants  C(f,‘Ue,  CjEe,  and  CEEo  are 
time-average  powers  (see  Section  6.3.1).  The  minus  sign  which  precedes  the 
E^i1e~  term  on  the  right-hand  side  of  (6.75)  makes  the  phase  of  C^,e  equal 
to  the  phase  of  the  coefficient  of  See  (6.35).  The  Ej‘xle~,  Ejfe~ , 
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and  Ejf°  terms  in  (6.75)  were  patterned  after  the  ,  and 

terms,  respectively. 


6.3.1  Time-Average  Power 

The  time-average  power  of  the  normalized  electric  field  Lz  /  \j  ZoiMeo 

of  (6.75)  is  given  by  a  formula  similar  to  (5.31).  This  power  consists  only 
of  the  time-average  powers  associated  with  the  E^\1e+ ,  ,  E \fe~,  and 

El\\°~  terms  on  the  right-hand  side  of  (6.75).  The  time-average  ^-directed 
power  associated  with  the  Eoife+  term  in  (6.75)  is  1  VV.  The  time-average 
— z-directed  powers  associated  with  the  Eo\  e“,  Eu~,  and  Eh  °~  terms  in 
(6.75)  are  IC™6]2,  | Cj^e\2,  and  IC^0)2,  respectively. 


6.3.2  The  Coefficient  C™e 

In  this  subsection,  C™e  of  (6.77)  is  expressed  in  a  form  suitable  for  calcula¬ 
tion.  In  (6.77),  S™e  is  given  by  (6.56)  in  which  the  quantities  zTiU2/c  and 
<j)b2  appear.  Comparing  (6.15)  with  (3.76),  we  see  that,  since  ar0i  < 


c 


01 


(6.84) 


where 


01 


is  G&  when  r  =  0  and  5  =  1.  Here,  6  is  either  TM  or  TE. 


Setting  r  —  0  in  (6.29),  we  obtain 


=  <e 


(6.S5) 


Substitution  of  (6.61),  (6.84),  and  (6.S5)  into  (6.56)  gives 


cTA/e 

‘-'01 


7=1 9=0 


TM 

<7 


01 


E  ^  ( jKE)- 

P+9^0 


(6.86) 


The  quantity  Z™e0  in  (6.77)  is  given  by  (6.69): 


77/9™ 

yTMco  _  VPOI 
Z01  -  u 


(6.87) 
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Substituting  (6.S6)  and  (6.87)  into  (6.77),  we  obtain 


C  =  -uEE  £  {c, 


(  yiTM  j (3™  L3 
TMe,-yTM  |  Vpq  e 


-y=l  9=0  p=0 

P+1^0 

(  V'yTEf>j0™L3 
pTMenTE  /  Vpq  e 
+  U01  ,P«  I  “ 


(6.88) 


where 


rTMe,iTM 

°01,w 


s~iTMe,"iTE 

°01,p9 


V  C  (^')^ 


(6.S9) 


(6.90) 


6.3.3  The  Coefficients  and  Cf^° 

In  this  section,  of  (6.78)  and  CjEo  of  (6.79)  are  expressed  in  forms 

suitable  for  calculation.  In  (6.78)  and  (6.79),  Sj \Ee  and  Sff0  are  given  by 
(6.58)  and  (6.59),  respectively.  The  quantities  zTE4/c  and  7 JEa  appear  in 
both  (6.58)  and  (6.59).  Comparing  (6.16)  with  (3.108),  we  see  that,  since 
x'n  <  ka, 

7T  “ ft*!  (6.91) 


where  [(S^  is  G^  when  r  =  1  and  s  =  1.  From  (3.58),  we  have 


lTnE  =  jfif- 


(6.92) 


Substitution  of  (6.61),  (6.62),  (6.84),  (6.91),  and  (6.92)  into  (6.5S)  and  (6.59) 
gives 


s™'  =  77=7  £  £  £  *,{  KE]„  #  (’V,™  +  (f ) 

~\]xn  —  1  7=19=0  P=o  v 

p+7^0 
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x  ( sin  <t>o\  igS4,l„<i^((f)'/;,™-^re)\ 

V  4>o  )  0\\a  * 


(6.93) 


cTEo  _  ~  1 

°11  —  7== 
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7==  EE  E  ^{[6H„^(«™  +  (f)v,TE) 

V^n  -  1  -,=1 9=0  P=o  v  v  0  7  y 

P+9#0 


(  sin<f>0\ 

1  . 

[Q?\ 

n  *U3 1 

((■ 

V  <f>0  j 

'  /? 

|T£. 

n  < 


} 


(6.94) 


The  admittance  YxlEeo  in  (6.78)  and  (6.79)  is  obtained  by  substituting 
(6.92)  into  [1,  eq.  (B.54)]: 


oTE 

•yTEeo  _  Pi  1 

11  kT) 


Substituting  (6.87),  (6.95),  and  (6.93)  into  (6.78),  we  obtain 

2  r  ( yiTM  j0^mL3 

QTEe  _  ^  -^2  ^  "I  cJEe’^™  V C 

_ i _ r»  _ rv 


7=1 9=0  p=0 

p+9#0 
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,TEe,iTM 
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An  ^/x^2  -  1 

pc  ( sinA\  Xu{G[4))n<t>bp4  \ 

+(*'V  4>o  )  ft Ea  J 
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2?r6 


cTEe.7T£  /£^ 


fit 


“7 


'sin  <ft0\ 

,  <^c  /  /?n£a 


} 


(6.95) 


(6.96) 


(6.97) 


(6.98) 
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Substituting  (6.87),  (6.95),  and  (6.94)  into  (6.79),  we  obtain 


s^TEo 

un 


EE  E  { 


~r=  l  <?= 0  p= 0 
P+9^0 


t ~iTEo,-yTM 
U11  ,P9 
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yjTMei(3™Lz 
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pc 
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\/xn  —  1 

sin^  a:'n(Gf^))ii^  j. 


<t>0 


0uE  a 


I2nb 
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/  tPC\  * 

<G*  >"*■ 


r 


+<T~j 

The  0p’s  in  (6.97),  (6.9S),  (6.100),  and  (6.101)  are  obtained 
=  1  in  (6.2S)-(6.31).  From  (2.8)  and  (2.9),  we  have 


so  that 


Assuming  that  <j)0  <  7r/2,  we  have 


b_ 
2  a 


(6.99) 

(6.100) 

(6.101) 

iy  setting 

(6.102) 

(6.103) 

(6.104) 
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Substitution  of  (6.103)  and  (6.104)  into  (6.28)— (6.31)  yields 


(6.105) 

(6.106) 

(6.107) 

(6.108) 
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Chapter  7 

The  Tangential  Electric  Field 
in  the  Apertures 


The  tangential  part  of  the  electric  field  in  the  left-hand  aperture  Ai  is  called 
z).  The  tangential  part  of  the  electric  field  in  the  right-hand  aperture 
Ai  is  called  z).  These  tangential  parts  are  given  by 

£l'h  W)  =  X  7=1,2  (7.1) 

where  up  is  the  unit  vector  in  the  p-direction  and  M_ ^  is  given  by  [1, 
eqs.  (2.11)  and  (2.12)].  Substituting  [1,  eqs.  (2.11)  and  (2.12)]  into  (7.1), 
we  obtain 

=  E  E  {i?,™ b .*0=)) 

p= 0  q—0 

P+Q^O 

+V£E  (u,  x  z))  }  (7.2) 

where  the  double  summation  is  truncated  as  in  (5.1).  In  (7.2), 

v;™  =  0,  p  =  0  or  q  =  0.  (7.3) 

Substitution  of  [1,  eqs.  (2.13)  and  (2.14)]  into  (7.2)  gives 

z)  =  ■u(-ir  (^) 
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p— 0  q— 0 
P+9^0 


+fe£  £  {v,;™^(!/7+,^)  +  V!JTEe™(^+,z+)}  (7.4) 


p=0  g=0 

P+9^0 


where  ?/')'+  and  2+  are  given  by  [1,  eqs.  (2.15)— (2.17)] 


Vy+  =  (2  -  7)^0  +  (-1  Y<f>x0  +  ^ 


1  ^ 

z+  =  z  +  - 


where,  from  (6.102), 


Now,  e™  and  e™  are  the  y-  and  ^-components  of  e™  given  by  (A. 10), 
and  and  ej^q  are  the  y-  and  2-components  of  eJ?E  given  by  (A. 23)  so  that 
(7.4)  becomes 

T7 '(Al),,  \  /  [b  fSm(l)o\ 


Er’faz)  = 


c  \  <f>0 


■5  5  </?(< 


p=0  q= 0 
P+9^0 


pny'r+  .  qxz+ 

cos  — 7 —  sin - 

0  c 


I  £n€(j  (  1 


77=0  0=0 

P+9it0 


.  pxy7+  <77r\r+ 

sin  — 7 —  cos - . 

b  c 


We  want  to  normalize  E\Ay\<p,z)  of  (7.S)  by  dividing  by  —  |ikoi‘Ue+|nns 
where  the  subscript  “rms”  denotes  the  root  mean  square  value  of  the  trans¬ 
verse  part  over  the  waveguide  cross  section  at  2  =  0.  Recall  that  £^’A/e+  is 
the  electric  field  of  the  2  traveling  wave  in  the  circular  waveguide.  We  choose 
to  divide  by  —  |£j1A/e+|rms  rather  than  |ikJAre+|rms  because  the  2-directed  elec¬ 
tric  current  associated  with  £jIA/e+  is  negative  at  (/>,  2)  =  (a,  0).  This  electric 
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current  is  given  by  —  up  x  Hj*{e+  where  H™e+  is  given  by  [1,  eqs.  (5.2)  and 
(5.9)]. 

According  to  (6.74),  the  transverse  part  of  Ec\Ie+  is  Z™eoe£ife{p,  4>) 
.e-i0oiM*  so  that 

i£jr,+ u  =  \z™‘°\ j [pdp  r  ■  £?"')}I/2  ■  <7-9) 

One  of  the  orthogonality  relationships  in  [1,  eq.  (B.64)]  and 

nB™ 

ZTMeo  _  7Poi_  (7.10) 

k 

reduce  (7.9)  to 

|£Sl/e+|rms  =  ('•!!) 

Equation  (7.10)  was  obtained  by  replacing  l/(u ;e)  by  rj/k  in  (6.69).  Substi¬ 
tution  of  (7.11)  and  (7.5)-(7.7)  into  (7.8)  gives 


3in  {  (^)  +  (2  “  ?)*)  +  y  }  cos  (~  +  t)  ’  (7-12) 


When  z  =  0,  the  ^-component  of  (7.12)  is  E^Al'\<^10)/j£LJ1Ee+jrms  given 
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by 


When  <j>  =  (2— 7)tt,  the  z-component  of  (7.12)  is  £(^((2— 7)71- ,  2)/|jE^/e+|rms 
given  by 


Chapter  8 

Numerical  Results 


A  computer  program  was  written  to  calculate  the  time-average  power  trans¬ 
mitted  into  the  rectangular  waveguides,  the  time-average  power  reflected  in 
the  circular  waveguide,  and  the  4>  and  2  components  of  the  electric  field  in 
the  apertures.  This  computer  program  will  be  described  and  listed  in  a  sub¬ 
sequent  report.  Some  numerical  results  obtained  by  means  of  this  computer 
program  are  presented  in  this  chapter. 

When  the  time-average  incident  power  in  the  circular  waveguide  is  unity, 
the  time-average  power  transmitted  into  the  rectangular  waveguides  is  Pt 
given  by  (5.49)  and  the  time-average  power  reflected  in  the  circular  waveguide 
is  called  PT.  According  to  the  discussion  in  Section  6.3.1, 

Pt  =  |C™'P  +  \C™r  +  |Cf,EY.  (8.1) 


Figures  8.1  to  8.5  show  plots  of  Pt  and  PT  versus  ka  for  various  values  of  L3 
when 


(8.2) 


The  plots  of  Figs.  8.1  to  8.5  are  for  L3/[\™]ka=2. 95  =  0.35,  0.40,  0.45,  0.50, 
and  0.55,  respectively.  Here,  [AoiMUa=2.95  is  the  wavelength  of  the  T M0l 
mode  in  the  circular  waveguide  when  ka  =  2.95: 


*a=2.95 


27TO 

^(2.95)2  -  xlx ' 


(8-3) 
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Fig.  8.1.  Plots  of  the  ratio  Pt  of  the  transmitted  power  to  the  input  power 
and  the  ratio  Pr  of  the  reflected  power  to  the  input  power  when 
L3  =  0.35  _2  .  The  input  power  is  the  power  of  the  inci¬ 

dent  TM0l  wave  in  the  circular  waveguide.  Pt  +  PT  =  1. 


Fig.  8.2.  Plots  of  the  ratio  Pt  of  the  transmitted  power  to  the  input  power 
and  the  ratio  Pr  of  the  reflected  power  to  the  input  power  when 
L3  =  0.40  (A^M]  .  The  input  power  is  the  power  of  the  inci¬ 

dent  TM0i  wave  in  tfie  circular  waveguide.  Pt  +  PT  =  1- 


S 


2.85  2.90  2.95  3.00  3.05  3.10 

Fig.  8.3.  Plots  of  the  ratio  Pt  of  the  transmitted  power  to  the  input  power 
and  the  ratio  PT  of  the  reflected  power  to  the  input  power  when 
L3  =  0.45  X™  .  The  input  power  is  the  power  of  the  inci- 
dent  T Mqi  wave  in  the  circular  waveguide.  Pt  -f  Pr  =  1. 


2.85  2.90  2.95  3.00  3.05  3.10 

Fig.  8.4.  Plots  of  the  ratio  Pt  of  the  transmitted  power  to  the  input  power 
and  the  ratio  PT  of  the  reflected  power  to  the  input  power  when 
L3  =  0.50  .  The  input  power  is  the  power  of  the  inci¬ 

dent  TM oi  wave  in  the  circular  waveguide.  Pt  +  PT  =  1. 
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Fig.  8.5.  Plots  of  the  ratio  Pt  of  the  transmitted  power  to  the  input  power 
and  the  ratio  Pr  of  the  reflected  power  to  the  input  power  when 
£3  =  0.55  [a™]  The  input  power  is  the  power  of  the  inci¬ 

dent  TMqx  w’ave  in  the  circular  waveguide.  Pt  +  PT  =  1. 


With  the  value  of  x01  given  in  [5,  page  2],  (8.3)  becomes 


ka=2.95 


2  na 

y(2^957^(2!404S2556)2 


3.67738806a. 


(8.4) 


The  value  ka  =  2.95  was  chosen  because  it  is  fairly  central  to  the  range  of 
values  of  ka  in  Figs.  8.1  to  8.5  (see  the  next  paragraph).  The  values  of  L\, 
and  L2  do  not  matter  because  the  loads  Z\  and  Z2  were  chosen  to  be  matched 
loads,  that  is, 

z.  =  z2  =  Z? 0E  (S.5) 

The  curves  of  Figs.  8.1  to  S.5  are  plotted  for  the  entire  range  of  ka  such 
that  only  the  TE]0  mode  propagates  in  the  rectangular  waveguides  and  only 
the  T A/01  and  T Eu  inodes  propagate  in  the  circular  waveguide.  Since  only 
the  TE io  mode  propagates  in  the  rectangular  waveguides, 


7r  <  kb  <  min 


(8.6) 
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where  “min”  denotes  the  minimum  of  the  values  in  parentheses.  Since  only 
the  T M0i  and  TEU  modes  propagate  in  the  circular  waveguide, 


x0i  <  ka  <  x2l.  (8.7) 

Substituting  (8.2)  and  (8.3)  into  (8.6),  we  obtain 

2.85599332  <  ka  <  5.7119866.  (8.8) 

Taking  the  values  of  z0i  and  x21  given  in  [5,  pages  2  and  32],  (8.7)  becomes 

2.40482556  <  ka  <  3.05423693  (8.9) 

Combining  (8.8)  and  (8.9),  we  have 

2.85599332  <  ka  <  3.05423693.  (S.10) 

In  Figs.  8.1  to  8.5,  Pt  approaches  zero  as  ka  approaches  2.S559932.  This 
is  expected  because  the  TEio  mode,  which  carries  the  transmitted  power, 
ceases  to  propagate  when  ka  becomes  less  than  2.S559932. 

The  numerical  data  of  Figs.  8.1  to  8.5  were  computed  with 

BKM  =  15  (8.11) 

XM  =  40.  (8.12) 


The  parameters  BKM  and  XM  are  not  written  with  italicized  letters  because 
they  are  input  variables  for  our  computer  program.  The  parameter  BKM  is 
introduced  in  (A.l)  and  used  in  (5.2).  The  constraint  (5.2)  on  the  values  of 
p  and  q  determines  the  order  of  the  moment  matrix  [T1  +  Y2  +  V'3]  which 
appears  in  [1,  eq.  (2.22)].  When  BKM  =  15,  the  order  of  the  moment  matrix 
is  32.  The  parameter  XM  is  introduced  in  (B.4).  The  effect  of  XM  is  to 
truncate  the  doubly  infinite  sum  that  appears  in  (3.1 )— (3.4).  The 

truncated  sum  is  I3r=ox  where  smax,  which  depends  on  r,  is  the  largest 

value  of  s  such  that 


jo,,  <  XM,  r  =  0 
j'\,,  —  XM,  r  =  1,2,-** 


(8.13) 


Assuming  that  XM  >  jo,i,  rmax  is  the  largest  value  of  r  such  that  j'T  l  <  XM. 
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Some  of  the  time-average  incident  power  is  transmitted  into  the  rectangu¬ 
lar  waveguides.  The  rest  of  it  is  reflected  in  the  circular  waveguide.  Therefore, 

Pt  +  Pr  =  1.  (8.14) 

The  plots  of  Pt  and  Pr  shown  in  Figs.  8.1  to  8.5  do  indeed  satisfy  (8.14). 
Because  the  ”alues  of  Pt  and  Pr  computed  separately  from  (5.49)  and  (8.1) 
satisfy  (S.12)  we  have  some  confidence  in  their  accuracy. 

The  magnetic  field  due  to  the  impressed  source  jZ‘mp  radiating  in  the 
circular  waveguide  with  the  apertures  closed  is  H}3\iLmp,Q.)  given  by  [1, 
eq.  (5.11)].  The  electric  current  at  z  =  0  on  the  wall  of  the  circular  waveguide 
associated  with  this  magnetic  field  is  J  given  by 

O  —jP^^Ls 

J  =  -uz  e  - cos(#fiiV/f3).  (S.15) 

\jKa 

The  magnitude  of  J  of  (8.15)  is  maximum  when  L3  =  0.5  A™  where  A™  is 
the  wavelength  of  the  TM01  mode  in  the  circular  waveguide.  If  an  aperture  is 
put  where  the  electric  current  would  otherwise  be  maximum,  the  tangential 
electric  field  will  be  large  in  this  aperture.  A  large  aperture  field  gives  a  large 
transmitted  power  Pt.  Thus,  we  expect  Pt  to  be  large  at  ka  =  2.95  in  Fig.  8.4 
because  L3  is  then  equal  to  0.5Ao  A/.  Actually,  Pt  is  even  larger  at  ka  =  2.95 
in  Figs.  8.2  and  8.3  where  L3  =  0.4A™  and  0.45AojA/,  respectively.  When 
ka  is  held  at  2.95,  the  curve  of  Pt  versus  L3/ shown  in  Fig.  8.6  attains  a 
maximum  at  a  value  of  L3/\™  somewhat  less  than  0.5.  This  phenomenon 
may  be  due  to  the  finite  extent  of  the  aperture  in  the  ^-direction.  Note  that 
Pt  =  0  at  L3/\™  =  0.25  in  Fig.  8.6.  This  is  expected  because  the  magnitude 
of  J  of  (8.15)  vanishes  when  L3/\™  =  0.25.  In  this  case,  the  aperture  has 
little  effect  because  there  is  no  flow  of  electric  current  to  stop  at  z  =  0. 

The  data  for  the  plot  of  Pt  versus  L3/\™  of  Fig.  8.6  were  computed  with 


(S.  16) 


- 

=  1.1 

a 

c 

=  0.5 
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ka 

=  2.95 

BI\M 

=  15 

XM 

=  40 

Zx 

=  Z2  = 

?TE 

^01 
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Fig.  8.6.  Plot  of  the  ratio  Pt  of  the  transmitted  power  to  the  input  power 
versus  L3/\™  when  ka  =  2.95.  The  input  power  is  the  power 
of  the  incident  T M0 1  wave  in  the  circular  waveguide. 


The  curve  in  Fig.  8.6  was  terminated  at  L3/ A™  =  1.5.  As  Z^/A™  becomes 
larger  and  larger,  Pt  versus  L3/\™  becomes  more  and  more  periodic  with 
period  one.  Adding  one  to  the  value  of  L3/\™  does  not  change  the  reflection 
of  the  TMqi  wave  from  the  short  at  z  —  L3.  It  only  changes  the  reflections  of 
the  even  and  odd  T E\\  modes  and  all  the  nonpropagating  modes.  Now,  with 
the  parameters  of  (8.16),  our  solution  for  the  electric  field  in  the  circular 
waveguide  did  not  contain  any  TE\\  modes.  The  computed  values  of  the 
constants  CjEt  and  Cjf°  in  (6.75)  were  zero.  When  L3/ is  large,  there 
is  not  much  reflection  of  nonpropagating  modes  from  the  short  at  r  =  L3 
because  any  nonpropagating  mode  suffers  attenuation  on  its  journey  from 
the  aperture  to  the  short  at  z  =  L3.  Furthermore,  its  reflection  suffers  the 
same  amount  of  attenuation  in  going  from  the  short  at  c  =  Z,3  back  to  the 
aperture. 

In  Fig.  S.6,  Pt  could  not  be  plotted  for  L3/ Ajju  <  0.067983034  because 
L3  cannot  be  less  than  c/2.  The  approximate  value  0.0679S3034  is  obtained 
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by  writing 


(8.17) 


and  substituting  (8.3)  and  (8.5)  into  the  right-hand  side  of  (8.17)  to  obtain 


L3  _  \/(2.95)2  -  (2.404S2556)2 

A0TiV/  St 


0.067983034. 


(8.18) 


Figures  8.7  to  8.13  show  plots  of  \E^A2\(f>,  0)\/\E^^Ie+  |rms  of  (7.13)  ver¬ 
sus  <t>/4>0  and  |£’('42)(0,  z)|/|F^’1A/e+|rrns  versus  z/c  for  ka  =  2.855994,  2.86, 
2.90,  2.95,  3.00,  3.05,  and  3.054236  when  Lzl[\™]ka-2$s  =  0.5.  As  given  by 
(7.13),  E^A2\<j>,  0) / |ikinMe+  Inns  is  the  ^-component  of  the  normalized  elec¬ 
tric  field  at  z  =  0  in  the  right-hand  aperture  A2.  As  given  by  (7.14), 
E(zA2){0,z)/\£SlMe+\rms  is  the  z-  component  of  the  normalized  electric  field 
at  (f>  =  0  in  A2.  The  values  2.855994  and  3.054236  were  purposely  chosen 
close  to  the  lower  and  upper  bounds  in  (8.10).  The  curves  in  Figs.  8.7  to  8.13 
are  not  smooth  because  they  were  obtained  by  drawing  straight  lines  between 
points  spaced  0.025  apart  in  4>/<i>o  and  0.05  apart  in  z/ zQ.  In  Figs.  8.7  to  8.13, 
\E[A2)(0iZ)\  is  generally  much  larger  than  |ik^A2*(<f>, 0)j.  This  is  expected  be¬ 
cause  the  aperture  A2  stops  only  the  r-directed  electric  current  J_  of  (8.15). 
There  is  no  ^-directed  electric  current  to  stop. 

The  data  for  the  plots  in  Figs.  8.7  to  8.13  were  computed  with 


a 


a 


L3  —  0.5[A™]*a=2.95  ' 
BI<M  =  33 


(8.19) 


XM  =  100 

Zt=Z2  =  Z lE 

where  [Ao{v/]fca=2.95  is  given  by  (8.5).  For  the  data  in  (S.19),  the  computed 
values  of  V/p))T'u  and  V^E  satisfy 


T/l  TX1  \/2TM 

pq  ¥pq 


i/i  TE  —  y2TE 
p< 1  pq 


(8.20) 
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The  ratios  (a)  |^2)|/|£0™e+Us  and  (b)  |£<">|/|££W*+U  of  the 
magnitudes  of  the  <j>-  and  z-directed  electric  fields  in  the  aperture 
A2  to  the  root  mean  square  value  of  the  electric  field  of  the  incident 
TMqi  wave  in  the  circular  cylinder  when  ka  =  2.S55994  and  I3  = 

°-5[«'L=,ss. 


Fig.  S.S.  The  ratios  (a)  |4",l/l£m"e+U  and  (b)  |£<'”>|/|fSw,+  U  of  the 
magnitudes  of  the  <f>-  and  ^-directed  electric  fields  in  the  aperture 
A2  to  the  root  mean  square  value  of  the  electric  field  of  the  incident 


TM01  wave  in  the  circular  cylinder  when  ka  =  2.86  and  L3  = 


1.95 


Fig.  8.9.  The  ratios  (a)  \E[*v\/\ £™e+|nn,  and  (b)  |£*”)|/l£&We+lr™  of  the 
magnitudes  of  the  d>-  and  z-directed  electric  fields  in  the  aperture 
T2  to  the  root  mean  square  value  of  the  electric  field  of  the  incident 
TMnX  wave  in  the  circular  cylinder  when  ka  =  2.90  and  L3  = 
0.5  |A™ 

01  Jfca=2.95 


Fig.  S.10.  The  ratios  (a)  l£i"’|/l£o"'+U  and  (b)  ''+i™.  of  the 

magnitudes  of  the  <t>-  and  z-directed  electric  fields  in  the  aperture 
A?  to  the  root  mean  square  value  of  the  electric  field  of  the  incident 
TiXIox  wave  in  the  circular  cylinder  when  ka  =  2.95  and  L3  = 


T 


1 

4>0 


Fig.  8.11.  The  ratios  (a)  \E^2)\/\^+U  and  (b)  |i^42)|/|££'fe+U  of  the 
magnitudes  of  the  <}>-  and  c-directed  electric  fields  in  the  aperture 
A2  to  the  root  mean  square  value  of  the  electric  field  of  the  incident 
the  circular  cylinder  when  ka  =  3.00  and  L3  = 


T M pi  wave  i 
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Fig.  S.12.  The  ratios  (a)  |£<-42,|/I£™‘+U  and  (b)  |£!J2»|/|£j,'''+|,™  of  the 
magnitudes  of  the  <f>-  and  r-directed  electric  fields  in  the  aperture 
A2  to  the  root  mean  square  value  of  the  electric  field  of  the  incident 
T A/01  wave  in  the  circular  cylinder  when  ka  =  3.05  and  L3  = 

0.5  [,\™l 


Jfca=2.95 


Fig.  8.13.  The  ratios  (a)  |4<!)|/I£™'+U  and  (b)  |£''1JI|/l£ori'''+U  of  tin 
magnitudes  of  the  d>-  and  c-directed  electric  fields  in  the  apertun 
A2  to  the  root  mean  square  value  of  the  electric  field  of  the  inciden 
T A/oi  wave  in  the  circular  cylinder  when  ka  =  3.054236  and  Z,3  = 

°"WL=,95- 


the  angle  7r  —  <^>  in  A2.  Furthermore,  the  z  component  of  the  electric  field 
at  the  angle  <j>  in  A\  is  equal  to  the  z  component  of  the  electric  field  at  the 
angle  n  —  (f>  in  A2. 

Note  that  the  values  of  BKM  and  XM  in  (8.19)  are  larger  than  those  in 
(8.16).  Larger  values  of  BKM  and  XM  were  needed  for  accurate  calculation  of 
the  tangential  electric  field  in  A 2  than  for  accurate  calculation  of  Pt.  Larger 
values  were  needed  because  the  convergence  of  the  electric  field  in  /42  with 
increasing  BKM  and  XM  was  slower  than  the  convergence  of  Pt.  The  electric 
field  in  A2  converges  slowly  because  its  ^-component  tends  toward  infinity  as 
1/(1  —  (4>I4>0)2)1'  when  4>/<f>0  approaches  ±1  where  u  =  (1 — 2</>0/ ?r)/ (3 — 2^0/ 7r), 
and  its  2-component  tends  toward  infinity  as  1/(1  —  (2 2/c)2)1/3  when  2 z/c 
approaches  ±1  [6,  page  387]. 
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Appendix  A 

Ordering  the  Expansion 
Functions 


The  expansion  function  A/^n(<ft.  z )  for  the  equivalent  magnetic  current  A/(I* 
is,  according  to  [1,  eq.  (2.13)],  defined  in  terms  of  esmn.  Here,  6  may  be  either 
TM  or  TE,  e™  is  the  electric-type  mode  vector  of  the  T Mmn  rectangular 
waveguide  mode,  and  is  the  electric-type  mode  vector  of  the  T  Emn  rect¬ 
angular  waveguide  mode.  The  expansion  function  M^n(6.  z)  is,  according  to 
[1,  eq.  (2.14)],  also  defined  in  terms  of  §?mn. 

The  expansion  functions  are  arranged  in  the  following  order: 

1.  {jWS",m=l,2,-.-,MM(2)},  {A*£",m  =  l,2,-,MM(3)}, 
1,2,.-.,  MM  (4)},-.. 

2-  {m2E,">  =  1,2,  — ,MM(1)}, 

{4*2®,  "•  =  0,1,-,  MM  (3)},- 

3.  {4*2", m  -1,2,-, MM (2)}  ,  {A£X",m  -  1,2,-,MM(3)}, 
{m2",  m  =  1, 2,  ■  •  • ,  MM  (4)}  ,  ■  -  • 

4.  {4fi',"i=l,2,-,MMil)},  {M2E."*=0,1,— ,MM(2)}, 

{m2e,™  =  o,i,---,mm(3)},---. 

As  listed  above,  items  3  and  4  are,  respectively,  items  1  and  2  with  the 
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superscript  1  replaced  by  2. 

The  Mil's  in  items  1,  2,  3,  and  4  are  variables  that  are  calculated  in  our 
computer  program.  The  value  of  MM(n  +  1)  is  determined  in  the  following 

{lT  1  r  lTKl  f  1 

ILmn  j  {Ji-mn  {-^mn 

and  {M2JE  }  are  terminated  by  requiring  that  the  nonnegative  integers  m 
and  n  be  small  enough  so  that 

kmnb  <  BKM  (A.l) 

where  kmn  is  the  mode  cutoff  wavenumber  given  by  (2.4),  and  BKM  is,  a 
specified  constant.  The  parameter  BKM  is  not  written  with  italicized  letters 
because  it  is  an  input  variable  for  our  computer  program.  Substitution  of 
(2.4)  into  (A.l)  gives 

<  BKM.  (A. 2) 

According  to  (A. 2),  MM  (n  +  1)  is  the  largest  integer  such  that 

<  BKM.  (A. 3) 

All  values  of  n  so  large  that  no  integer  MM  (n  4-  1)  satisfies  (A. 3)  are  disal¬ 
lowed. 

If  we  define  M.}TXf  to  be  the  ALlTXt  in  item  1  for  i  =  1,2,  •  •  •,  then 
there  will  be  a  pair  of  integers  m  and  n  such  that  M~  ‘  —  AL^n  ■  Thus, 

there  will,  for  the  be  a  correspondence  between  eacli  pair  of  integers 

m  and  n  in  use  and  eacti  single  integer  i.  Replacing  m,  n.  and  i  by  p.  q ,  and 
j,  respectively,  t^ere  will,  for  the  4/iru,s.  be  a  correspondence  between  each 
pair  of  integers  p  and  q  in  use  and  each  single  integer  j.  This  is  what  is 
called  in  the  sentence  following  (2.5)  the  correspondence  between  each  pair 
of  integers  (p,  q)  used  in  [1.  eq.  (3.44)]  and  the  subscript  j  in  [1.  eq.  (3.44)]. 

If  we  define  M},TE  to  1  e  the  AJ_lTE  in  item  2  for  i  =  1.2.  •  •  • .  then 
there  will  be  a  pair  of  integers  m  and  n  such  that  \f_)TE  =  Thus, 

there  will,  for  the  Al}TE' s,  be  a  correspondence  between  each  pair  of  integers 
m  and  n  in  use  and  each  single  integer  i.  Replacing  m.  n,  and  i  by  p.  q. 
and  j,  respectively,  there  will,  for  the  4£lTE  s,  be  a  correspondence  between 
each  pair  of  integers  p  ?nd  q  in  use  and  each  single  integer  j.  This  is  what  is 
called  in  the  sentence  following  (2.5)  the  correspondence  between  each  pair 
of  integers  (p,  q)  used  in  [1.  eq.  (3.47)]  and  the  subscript  J  in  [1,  eq.  (3.47)]. 


76 


Appendix  B 

Roots  of  Bessel  Functions  and 
Their  Derivatives 


This  appendix  describes  how  numerical  values  of  {jns,  s  =  1,2,  ••  •  ,5max;  n  = 
0, 1,  •  •  •}  and  {j’ns,  s  =  1,2,  •  •  •  ,smax;  n  =  0, 1,  •  •  •}  are  obtained.  Here,  jn,  is 
the  s1*1  root  of  Jn: 

Jn(jn,)=0.  (B.l) 

Furthermore,  j'n.  is  the  root  of  J': 

J'ntin,)  =  °-  (B-2) 


In  (B.l),  Jn  is  the  Bessel  function  of  the  first  kind  of  order  n.  In  (B.2),  J'n  is 
the  der:  'ative  of  Jn  with  respect  to  its  argument.  Now,  the  roots  {jns }  and 
{j'n,}  are  ordered  so  that 

0  <  jo,i  <  j'o  i  <  jo,2  <  jo, 2  ’  •  ’  <  jo  ,$m*x  <  Jo 

n  ^  in,  1  <  in,l  <  in, 2  ^  in, 2  <  in,3m»i  <  in.J i  U  =  1, 2,  •  ■  • 

Here,  smax  depends  on  n.  Given  n,  smax  is  the  largest  value  of  s  such  that 


io.3  <  XM,  n  =  0 
in,,  <  XM,  n  =  1,2,--* 


(B.4) 


where  XM  is  not  written  with  italicized  letters  because  it  is  an  input  variable 
for  our  computer  program.  Assuming  that  XM  >  j0,i,  aH  values  of  n  so  large 
that  j'n  j  >  XM  are  disallowed.  In  (B.  l),  XM  is  a  constant  that  controls  the 


number  of  roots  that  are  calculated.  If  one  wants  to  calculate  more  roots, 
one  should  choose  a  larger  value  of  XM. 

Our  numerical  values  of  s  =  1,2,  •••,49;  n  =  0,1,  •••,19}  and 
{j'ns  ,s  =  1,2,  •••,49;  n  =  0,1,  ••*,19}  are  taken  directly  from  [5,  Tables 
1  and  2].  For  other  values  of  n  and  s ,  we  calculate  jns  and  j'ns  by  means  of 
formulas  given  in  [5].  In  the  body  of  the  present  report,  the  sth  roots  of  Jn 
and  J'n  were  called  xns  and  x'ns  to  coincide  with  the  notation  in  [1].  H<jre  in 
Appendix  B,  these  roots  are,  more  in  harmony  with  the  notation  in  [5],  called 
jns  and  j'ns.  However,  our  notation  in  Appendix  B  is  slightly  different  from 
that  of  [5].  Our  jns  is  what  Olver  calls  j „)S.  Our  j'0s  is  what  Olver  calls  jg  +1. 
Our  j'ns  for  n  >  1  is  what  Olver  calls  j'n  s.  Our  “inner  fringe”  calculated  values 
of  {i2os,  s  =  1,2,  •••,50},  {j„5o,  n  =  6, 1,2,  •••,19},  {j'20s,  s  =  1,2,  -  -  -  ,50}, 
and  {inso»  n  =  0, 1,2,  •  •  • ,  19}  agree  well  with  the  “outer  fringe”  tabulated 
values  in  [5,  Tables  1  and  2].  There  is  no  entry  in  [5,  Table  2]  which  corre¬ 
sponds  to  our  j'0  50.  However,  according  to  [3,  formula  9.1.2S],  j'0  50  =  jli50, 
and  there  is  an  entry  in  [5,  Table  1]  which  corresponds  to  our  j150.  Although, 
as  stated  earlier  in  this  paragraph,  our  notation  does  not  normally  place  a 
comma  between  the  subscripts  of  j  and  j we  had  to  use  a  comma  in  the 
previous  two  sentences  in  order  to  separate  the  “0”  from  the  “50”  and  the 
“1”  from  the  “50”. 


B.l  Evaluation  of  Roots  of  Bessel  Functions 
of  Large  Order 


For  n  >  20  and  s  >  1,  we  approximate  jn,  by  [5,  eq.  (9.01)] 


•  _  ,  Pi  .  P2 

Jns  —  T  T 

n 


where  z,  pi ,  and  p2  are  tabulated  functions  of  —  (  where 

-  C  = 

In  (B.6),  a,  is  the  negative  root  of  the  Airy  function  Ai: 


Here, 


Ai(fl3)  =  0,  s  =  1,2,  ••  • . 

0  <  —a i  <  — «2  <  —03  <  •  • 


(B.5) 


(B.6) 


(B.7) 

(B.8) 


78 


The  roots  {a3 ,  s  =  1,2,  •••,50}  are  tabulated  [5,  Table  V,  page  78].  For 
5  >  50  [5,  eq.  (9.07)], 

“■=-A2,3(i+4*y  (m) 


where  [5,  eq.  (9.09)], 


A  =  |(4S-1). 


(B-10) 


Actually,  z,  px,  and  P2  are  tabulated  functions  of  — (  only  for  (0  <  —  £  <  7.5) 
[5,  Table  IV,  pages  72  and  74].  If  —  £  >  7.5,  then  z  —  2/(3£3)  and  p1  are 
tabulated  functions  of  £  [5,  Table  IV,  page  74]  where 


£  = 


1 


V=(' 


(B.li) 


If  —(  >  7.5,  then  p2  —  0. 

The  modified  interpolation  formula  of  Everett  [5,  eq.  (9.04)],  [7,  page  57], 


/p  =  (1  -  p)fo  +  pfi  +  E262m0  +  +  Mj7o  +  JV47?, 


(B.12) 


is  used  to  obtain  accurate  values  of  2,  pi,  and  p2.  In  (B.12),  the  f's  are  vaiues 
of  the  function  being  interpolated,  the  6^’s  are  modified  second-order  differ¬ 
ences,  and  the  74,s  are  modified  fourth-order  differences.  In  particular,  fp  is 
the  interpolated  value  of  /  at  the  actual  value  of  the  argument;  /0,  and 
7o  are  the  tabulated  values  of  /,  8^,  and  74  at  the  nearest  smaller  tabulated 
value  of  the  argument;  and  /l5  and  74  are  the  tabulated  values  of  /,  6^, 
and  74  at  the  nearest  larger  tabulated  value  of  the  argument.  The  argument 
is  either  —(  or  £.  Numerically,  p  is  the  ratio  of  the  difference  between  the 
actual  value  and  the  nearest  smaller  tabulated  value  of  the  argument  to  the 
difference  between  the  nearest  larger  and  nearest  smaller  tabulated  values  of 
the  argument.  Thus,  the  actual  value  of  the  argument  is  a  fraction  p  of  the 
way  from  the  nearest  lower  tabulated  value  to  the  nearest  upper  tabulated 
value.  In  (B.12),  F2,  F2 ,  A/4,  and  N,x  are  given  by  [7,  pages  56  and  57] 


E2 

f2 

AU 


p(p  ~  l)(p  ~  2) 

6 

(p  +  l)p(p  -  1) 

6 

1000  E2  I <P  +  l^~  —  +  0.184 


(B.1.3) 

(B.14) 

(B.15) 
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N4  =  lOOOi^  |  —  2^P . -2)-  +  0.184| .  (B.16) 


B.2  Evaluation  of  Roots  of  Derivatives  of 
Bessel  Functions  of  Large  Order 


For  n  >  20  and  s  >  1,  we  approximate  j'ns  by  [5,  eq.  (9.02)] 


*/ 

•7  ns 


.  ?i  .  ?2  qz 

nz  + - h  —z  +  -r 

n  n3 


(B.17) 


where  z,  q4,  q2 ,  and  <?3  are  tabulated  functions  of  —  £  [5,  Table  IV,  pages 
72-75]  where 

-C  =  -n-2/V,.  (B.18) 

In  (B.18),  a,  is  the  s^1  negative  root  of  Ai\  the  derivative  of  the  Airy  function 
Ai: 

AiK)  =  0,  5  =  1,2,-..  .  (B.19) 


Here, 


0  <  -a[  <  -a2  <  -a'3  <  •  •  • . 


(B.20) 


The  roots  {a',  s  =  1,2,  •••,50}  are  tabulated  [5,  Table  V,  page  78].  For 
s  >  50  [5,  eq.  (9.08)], 

<  =  -"2/3  (‘  -  w)  (B'2,) 

where  [5,  eq.  (9.09)] 

M  =  y(4s-3).  (B.22) 

The  tabulation  of  2  was  described  in  the  two  sentences  following  (B.10). 
Actually,  q4,  q% ,  and  <73  are  tabulated  functions  of  —  £  only  for  0  <  £  <  7.5 
[5,  Table  IV,  pages  73  and  75].  If  —  £  >  7.5,  then  q2  =  <73  =  0  and  qi  is 
a  tabulated  function  of  £  [5,  Table  IV,  page  75]  where  £  is  given  by  (B.ll). 
The  interpolation  formula  (B.12)  is  used  to  obtain  accurate  values  of  2,  qx, 
<72,  and  q3. 


SO 


B.3  Evaluation  of  Large  Roots  of  Bessel 
Functions 

For  0  <  n  <  19  and  s  >  50,  we  approximate  jna  and  j'ns  by  the  truncated 
McMahon  expansions  [5,  eqs.  (1.10)  and  (1.12)] 

4 

__  a  V''  "T-2r-l  /  r, 


b-Y - - 

(2r  -  l).'23'/'"1 


j'ns  =  &  -  S 


/i2r-l 


In  (B.23), 


(2r  -  l)!23^'<2r-i) 


0  =  (2n  +  4s  -  1)- 
4 


and  [5,  eq.  (1.11)] 

Ax  =  fi  -  1 

A3  =  (//  -  l)(7/t  -  31) 

A5  =  4(/i  -  1)(83/  -  982/i  +  3779) 


(B.23) 

(B-24) 

(B.25) 

(B.26) 

(B.27) 

(B.28) 


At  =  6 (ft  -  l)(6949/t3  -  153855/  +  1585743 p  -  6277237)  (B.29) 


where 


/x  =  4  n2.  (B.30) 

In  (B.24), 

/  =  (2n  +  4s-3)j  (B.31) 

and  [5,  eq.  (1.13)] 

A[  =  ft  +  3  (B.32) 

4  =  7/z2  +  82/x  -  9  (B.33) 

A's  =  4(S3/  +  2075/x2  -  3039 p  +  3537)  (B.34) 

A'7  =  6(6949/  +  296492/  -  1248002/  +  7414380/x  -  5853627)  (B.35) 
where  /i  is  given  by  (B.30).  Expressions  (B.23)  and  (B.24)  expand  to 


]n3  ~  13  8 (3  384/  15360/  (  4  )  3440640/  (  6  )  *  '  ^ 

,  A\  A'3  1  M'\  1  //\ 

Jna  ~  J  80'  384/3  15360/?' 5  V  4  )  3440640/7  l  6  )  '  ( 
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Appendix  C 

The  Effect  of  Loads  on  the 
TMqi  Source  in  the  Circular 
Waveguide 


In  the  body  of  this  report,  we  find  the  electric  fields  in  the  waveguide  mode 
converter  due  to  the  excitation  of  a  unit-amplitude  2-traveling  T MQX  wave. 
The  source  of  this  wave  is  taken  to  be  the  electric  current  source  J_'mp  whose 
-2-traveling  waves  see  a  matched  load.  Specializing  further,  we  take  J_'mp 
to  be  the  transverse  electric  surface  current  density  located  at  2  =  —  L5  that 
launches  the  2-traveling  wave  whose  electromagnetic  field  is  ( Ej^fe+ .  ) 

in  the  region  for  which  2  >  —  Z5.  Here, 


rpTMe- 1- 
£±0l 


=  4>)e~^  2  +  u 


£01 


TMe+ 


01 
TMc 


^01 


(p,<t>)e  30o 


TM  , 
1  ' 


jut 


.  (C.l) 


Equation  (C.l)  was  obtained  by  substituting  (6.63)  into  [1,  eq.  (B.  1 )].  The 
electric  and  magnetic  fields  defined  by  (C.l)  are  those  of  the  2-traveling 
T A/01  mode.  We  assume  that  J_>mp  also  launches  the  -2-traveling  wave  whose 
electromagnetic  field  is  C(E%xfe~,  H^XU~)  in  the  region  for  which  2  <  —L$. 
Here,  C  is  an  unknown  constant  and 


TM.  x 


FT\fe-  _ 
£^01  ~ 

riTMe-  _ 

U-01  ~ 


(p,<f>)eJ0™2 


TM  _ 


+ 


{^mf)2^oi,e(p^y0o{  z 


1  T  Me 
£01 


jut 


■  (0.2) 
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Equation  (C.2)  was  obtained  by  substituting  (6.63)  into  [l,  eq.  (B.2)].  The 
electric  and  magnetic  fields  defined  by  (C.2)  are  those  of  the  -2-traveling 
TMq\  mode.  Requiring  the  transverse  electric  field  to  be  continuous  at  2  = 
—  L$,  we  obtain 

C  =  -e*0™L\  (C.3) 


The  TA/0i  waves  radiated  by  J‘mp  in  the  circular  waveguide  are  shown  in 
Fig.  C.l  where  “1”  is  the  coefficient  of  the  2-traveling  TA/0i  mode  in  the 
region  for  which  —L$  <  z  <  —c/2.  ‘T”  is  the  coefficient  of  the  —2-traveling 
TA/01  mode  in  the  region  for  which  —  L$<z<  —c/2.  “C  +  T”  is  the 
coefficient  of  the  —  ^-traveling  TA/01  mode  in  the  region  for  which  —  co  <  z  < 
—  T5.  The  TA/01  waves  shown  in  Fig.  C.l  are  those  dealt  with  in  the  body  of 
this  report.  Thus,  from  (6.75),  we  have 


T  =  -Clxxiee~l20™L:i. 


(C.4) 


Here,  T  is  a  reflection  coefficient  because  7  is  the  ratio  of  the  coefficient  of  the 
-^-traveling  mode  field  (C.2)  to  the  coefficient  of  the  2-traveling  mode  field 
(C.l)  in  the  region  for  which  —  L5  <  2  <  —c/2.  Since  the  magnetic  fields  of 
the  mode  fields  are  h^xte(p,  rather  than  ±hTAfe(p,  <p)e*]i3™z,  7  is 

a  reflection  coefficient  for  the  current  rather  than  for  the  voltage. 

The  loads  mentioned  in  the  title  of  this  appendix  are  taken  to  be  the 
TA/01  loads  Zi 4  at  2  =  —LA  and  Zi6  at  2  =  —  L6.  See  Fig.  C.2.  A  TA/01  load 
is  a  load  that  acts  on  the  voltage  and  current  of  only  the  T A/0 1  waves.  When 
the  loads  Zn  and  Zl§  are  in  place,  the  TA/0i  electric  and  magnetic  fields  K 
and  //  in  the  circular  waveguide  are,  as  indicated  in  Fig.  C.2,  assumed  to  be 
given  by 


r*  _  /°»+  rpTMe- f-  >  /^—  rpTMe— 

Ll  —  ^4  -r  U4  £j_qX 

K  =  CillZ?u+  +  C;hZ?'- 

p  —  p+  pTMe+  1  p>-  pTMe- 
El  —  o5  £jqX  -r  o5  £lqX 

K  =  CtHiP*  +  C;hZ?"- 


P  —  pTMe+  I  c-rTA/e- 
El  —  ^6  £4)1  r  C-g  £lqx 

&  =  (%#£*•*  +  Ca-fg il'~ 


<  2-  <  L 5 


—  L  5  <  2  <  —L$ 


-  U<z  <  -- 


(C.5) 


(C.6) 
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«-  c  +  r 


1  -> 


r 


+—  to  z  =  —  oo 


l3 


Fig.  C.l.  The  T Mqi  waves  radiated  by  Jimp  in  the  circular  waveguide. 

The  situation  in  Fig.  C.l  is  the  same  as  that  in  the  body  of 
this  report.  There  are  no  additional  loads. 


Fig.  C.2.  The  T  Mqi  waves  radiated  by  j/imp  when  the  additional  loads 
ZL4  and  Zl6  are  present  in  the  circular  waveguide. 


S4 


where  the  C's  are  unknown  constants. 

The  effect  of  the  loads  Z^4  and  is  to  change  the  electric  field  inci¬ 
dent  on  the  aperture-perforated  section  of  circular  waveguide  from  1'We+ 
in  Fig.  C.l  to  Cq  MoiIe+  in  Fig.  C.2.  Viewing  this  incident  electric  field  as 
the  excitation  of  the  aperture-perforated  section  of  circular  waveguide,  we 
deduce  that  the  fields  in  this  section  of  waveguide  in  Fig.  C.2  are  those  in 
Fig.  C.l  multiplied  by  C$.  Therefore, 

ci  =  r  c+.  (c.8) 

Equation  (C.8)  is  one  simultaneous  equation  in  the  variables  C4,  C4,  C+, 
C$  ,  Cg,  and  C6".  In  the  following  five  paragraphs,  we  obtain  five  more 
simultaneous  equations  in  these  six  variables. 

Because  the  waves  launched  by  JItnp  in  Fig.  C.2  are  the  same  as  those  in 
Fig.  C.l,  we  have 


Ct  =  Ct  +  1  (C.9) 

c4  =  C5*  +  c.  (C.10) 


Recall  that  the  loads  Zn  and  Zlg  act  on  the  voltages  and  currents  of  the 
TMoi  waves.  These  voltages  and  currents  are  called  the  TMoi  voltages  and 
currents.  Seeking  to  define  the  TMoi  voltages  and  currents,  we  substitute 
(C.l)  and  (C.2)  into  (C.5)-(C.7)  and  take  only  the  transverse  part  of  E_, 
which  is  called  £(: 


K,  =  4,) 

H  =  {cte~’0™‘  +  C4-ei«“')  hln'ip-  <t>) 

E,  =  -  Cfe*®"') 

K  =  (< Cte ->«£"»  +  C5-erf*'=)  gf'ip,  4) 

£■  =  [C*e~>W-  - 

H  =  (C6+e-'<"'  +  C; *#•“’)  i 


-  L4  <  z  <  -L5  (C.ll) 


-  Ls  <  z  <  -Lg  (C.12) 


-L6<z<-C~.  (C.13) 


The  TMoi  voltages  are  defined  to  be  the  coefficients  of  e£fIe(p,<t>)  in  (C.ll)- 
(C.13).  The  TMqi  currents  are  defined  to  be  the  coefficients  of  h^^fe(p,  tj>)  in 
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(C.11)-(C.13).  Viewing  Fig.  C.2  as  a  circuit,  the  TM0 1  voltage  is  the  voltage 
of  the  upper  line  with  respect  to  the  lower  one,  and  the  TMq\  current  is  the 
z-directed  electric  current  on  the  upper  line. 

The  presence  of  Zi 4  at  z  =  —  L4  requires  that 

yl  .  =  -Zl  4  (0.14) 

where  V  is  the  TM0 1  voltage  and  I  is  the  TMq\  current.  The  subscript 
“z  =  —  L4”  denotes  the  limit  as  z  approaches  —  L\  from  above.  “From 
above”  means  through  values  which  are  greater  than  —  L4.  Extracting  V  and 
I  from  (C.ll)  and  substituting  them  into  (C.14),  we  obtain 


C4V'4  -  C;e-Jl'  ZL4 

C4  ejl*  +  C4~e~^  ~  Z0™“> 

(C.15) 

where 

i,  =  13™  L„ 

Solving  (C.15)  for  C4  in  terms  of  C4  ,  we  arrive  at 

(C.16) 

ct  =  r4c4- 

(C.17) 

where 

yTMeo  y 

r  _  ^01  z'L4  j2ii 

l*-z™*°  +  zL4e  ■ 

Since  there  is  no  series  load  at  z  =  —  Z,6, 

(C.18) 

(V'U-t.-  = 

(C.19) 

where  the  subscript  “z  =  —  L$  ”  denotes  the  the  limit  as  z  approaches  —  Z6 
from  below  and  “z  =  —  Lq"  denotes  the  limit  as  z  approaches  —  L6  from 
above.  Extracting  the  V’s  from  (C.12)  and  (C.13)  and  substituting  them 
info  (C.19),  we  obtain 

C5V'8  -  C5-e"J'6  =  C6V'6  -  C; e~J'6  (C.20) 

where 

le  =  (C.21) 
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At  2  =  -L6,  the  TMqi  current  V/ZL&  flows  from  the  upper  terminal  of  ZL6 
to  its  lower  terminal  so  that 


VU-L:  =  -^lvU-L;  +  li]„.Lr 


(C.22) 


Extracting  V  and  /  from  (C.12)  and  (C.13)  and  substituting  them  into 
(C.22),  we  obtain 


C,V7#  +  C7e~>k  =  Ctejl«  +  CZe~]k  + 


Cfe’1*  -  C5‘e-J'6 


(C.23) 


where 


Equation  (C.23)  becomes 


yT Mco  9 
^01 


(C.24) 


C5+(l  -  Z)eJ,«  -  Cf(l  +  Z)e~]l 6  +  <C6+ZeJ'6  +  =  0.  (C.25) 

Equations  (C.17),  (C.10),  (C.9),  (C.20),  (C.25),  and  (C.S),  ordered  as 
cited,  are  written  in  matrix  form  as 


i  -r4  0 

0  1  0 

-1  0  1 

0  0  —  eji® 

0  0  (l-Z)e-''6 

0  0  0 


-1  0 

0  0 


o  o  o  ire; 

-i  o  o  c ; 

0  0  0  Ct 

e~jl*  ejl6  —e~:l6  C$ 

-(1  +  Z)e~]l*  Zejlf>  Ze~]l*  C6+ 

o  -r  i  cz 


(C.26) 


We  proceed  to  solve  the  matrix  equation  (C.26)  for  C+ ,  C4  ,  C$ ,  C5  ,  C£ , 
An  d  Cg  . 
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Adding  the  first  row  of  (C.26)  and  the  product  of  T4  with  the  second  row 
to  the  third  row,  we  obtain,  in  view  of  (C.3), 


■  1 

-r4 

0 

0 

0 

0  ' 

\ct  1 

0 

i 

0 

-1 

0 

0 

c4- 

0 

0 

1 

-r4 

0 

0 

c5+ 

0 

0 

_  ej4 

e-]U 

e3k 

— e-jie 

c5- 

0 

0 

(1  -  Z)ejl 6 

-(1  +  Z)e~jt* 

Ze3' 6 

Ze~3k 

C6+ 

.  0 

0 

0 

0 

-r 

1 

.  ^6  . 

'  0  ' 
C 
D 
0 
0 
0 


(C.27) 


where 

d  =  i  -  r4ej2p™Li. 
The  last  four  rows  of  (C.27)  are 


(C.28) 


1 

— 

0 

0  ‘ 

\c?  1 

’  D  ' 

—e3k 

e3'6 

-e~jk 

C5' 

0 

(1  -  Z)e3‘* 

-(1  +  Z)e~3k 

Ze3‘« 

Ze~3k 

c6+ 

0 

0 

0 

-r 

1 

lc;  \ 

.  0  . 

(C.29) 


Adding  the  product  of  e]l<i  with  the  first  row  to  the  second  row  of  (C.29)  and 
adding  the  product  of  —(1  —  Z)e3‘6  with  the  first  row  to  the  third  row,  we 
obtain 


‘  1 

-r4 

0 

0  ' 

\ct 1 

0 

e~3le  -  r4eJ'6 

e3k 

-e~jk 

c5- 

0 

-(1  +  Z)e~3h  +  (1  -  Z) r4e* 

Ze3‘*  Ze~3k 

Q+ 

.  0 

0 

-r 

1 

Lq- 

D 

De3‘6 

-(1  -  Z)DeJl* 

0 


(C.30) 
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Discarding  the  first  row  of  (C.30)  and  multiplying  the  second  and  third  rows 
by  e~j!e ,  we  obtain 


e-j2'6  -  r4 

1 

—  e  i2le 

\C;  1 

-(l  +  Z)e->21°  +  (l  -  Z)T4 

z 

Ze~j2‘* 

C6+ 

0 

-r 

1 

.  ^6  . 

D 

-(1  ~Z)D 
0 


(C.31) 


Adding  the  product  of  the  first  row  with  {(1  +  Z)e  }2lt  —  —  Z)Y  4]  I  {e  j2/6  — 

r4}  to  the  second  row  of  (C.31),  we  obtain 


'  e~]2h  -  r4 

1 

-e->21* 

'c8-  • 

0 

(1+2 Z)e~i2l«  -  F4 

-e-j4'6  -  (2 Z  -  l)r4e^2'6 

C6+ 

e~j21*  -  r4 

e-J21*  -  r4 

0 

-r 

1 

.C6- . 

r  ^ 

1 

2  Z£>e~j2'6 

e -  r4 

0 


(C.32) 


Discarding  the  first  row  of  (C.32)  and  multiplying  the  second  row  by  e  j2le 
r4,  we  obtain 


■  (1  +  2Z)e-j2i«  -  r4 

_e-^'6  -  (2 Z  -  l)r4e-j2'6  ' 

■c6+- 

-r 

1 

.  ^6  . 

‘  2 ZDe~j21* 
0 


(C.33) 


Adding  the  product  of  the  first  row  of  (C.33)  with  T/{(1  +  2 Z)e~]2li  —  T4} 
to  the  second  row,  we  obtain 


1  - 


r(e"*'«  +  (2 Z  -  l)r4e--'2,«) 

(1  +  2Z)e~#‘*  -  r4 


cz  = 


2ZDYe~ j2'6 


(1  +2 Z)e~^  -  r4 


(C.34) 
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Solving  (C.34)  for  C6  ,  we  obtain 


2ZDT 
6  “  A 

where 

a  =  i  +  2Z  -  (2Z  -  i)rr4  -  r4ej2'6  -  re-j2'6. 

Substitution  of  (C.35)  into  the  second  row  of  (C.33)  gives 


(C.35) 

(C.36) 


(C.37) 


Substituting  (C.35)  and  (C.36)  into  the  first  row  of  (C.32)  and  solving  for 
C$  ,  we  arrive  at 


{(2Z-l)T  +  e*l*}D 
—  - 


(C.3S) 


Substitution  of  (C.3S)  into  the  first  row  of  (C.30)  gives 


{2Z  +  l-Ve-*l*}D 

— 


(C.39) 


Next,  C4  is  given  by  (C.10).  Finally,  C4  is  given  by  (C.17).  It  can  be  verified 
that  the  C’s  given  by  (C.35),  (C.37),  (C.3S),  (C.39),  (C.10),  and  (C.17)  do 
indeed  satisfy  (C.26). 
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